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1.  INTRODUCTION 


When  a  two-dimensional  boundary  layer  encounters  an  obstacle  protruding 
from  the  floor,  the  flow  experiences  the  change  of  the  pressure  field  near  the  obstacle. 
The  blockage  effect  due  to  the  obstacle  causes  the  approaching  mean  boundary  layer, 
initially  with  only  spanwise  vorticity,  to  be  deflected  to  also  have  a  longitudinal 
component  of  vorticity  (secondary  flow  of  Prandtl’s  first  kind).  Upstream  of  the 
obstacle,  the  flow  separates  from  the  floor  due  to  the  adverse  pressure  gradient.  The 
flow  entrained  down  along  the  edge  of  the  obstacle  is  recirculated  into  the  separated 
flow  region.  The  separation  with  the  secondary  flow  forms  a  very  complicated, 
three-dimensional  vortical  structure  in  this  region.  This  vortical  structure  trails 
downstream  along  each  side  of  the  obstacle.  A  main  feature  of  this  vortex  system 
is  its  persistence  downstream  (Sedney  and  Kitchens,  1975).  The  vortex  system  is  a 
form  of  “horseshoe”  vortex. 

The  formation  of  the  horseshoe  vortex  system  is  an  inviscid  process,  not  pro¬ 
duced  directly  by  viscosity  or  turbulent  stresses.  It  is  formed  around  any  type  of 
junction  whether  the  flow  is  laminar  or  turbulent  (Baker,  1978;  1979).  Yet,  in  a  tur¬ 
bulent  boundary  layer,  the  turbulent  Reynolds  stresses  also  create  the  cross-stream 
gradients  in  the  flow,  leading  to  the  secondary  flow  (Prandtl’s  second  kind).  Thus, 
the  flow  behavior  becomes  more  complicated  with  the  time-dependent  nature  of  a 
turbulent  flow. 

The  junction  vortex  flow  can  be  found  in  many  practical  examples,  such  as  the 
wing-fuselage  and  wing-pylon  junctions  of  aircraft,  the  appendage-hull  jimction 
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of  ships,  and  blade-shroud  junction  of  turbomachinery.  This  flow  structure  has 
undesirable  effects  on  the  flow  characteristics  in  many  cases;  degraded  aerodynamic 
performance,  increase  in  heat  transfer  in  this  region,  noise,  and  interference  with 
devices  that  follow  the  wing. 

Drag  increases  in  the  region  of  the  separated  flow  due  to  the  increased  shear 
stresses  caused  by  the  high  momentum  fluid  entrained  from  the  outer  part  of  the 
boundary  layer  flow  or  from  the  potential  flow  down  along  the  edge  of  the  obstacle. 
The  flow  around  a  bridge  pier  in  a  river  is  an  example  that  shows  the  increase  of  the 
surface  shear  stresses.  This  causes  accelerated  erosion  of  the  pier  base  as  the  sand 
upstream  of  the  pier  is  dug  out  and  flows  downstream.  Non-uniformity  of  the  flow 
in  the  wake  region  (Hawthorne,  1954)  interferes  with  the  operation  of  the  control 
devices  mounted  downstream.  In  the  junction  between  the  sail  appendage  and  the 
hull  of  a  submarine,  this  vortex  results  in  noise  due  to  the  time-dependent  nature 
of  the  vortical  flow  trailing  downstream. 

Attempts  have  been  made  to  reduce  these  undesirable  effects  caused  by  the 
junction  vortex.  Leading  edge  strakes,  and  fairings  or  fillets  have  been  used  to 
improve  the  aerodynamic  performance  (Bertelrud,  1988;  Devenport  ei  o/.,  1989). 
However,  to  achieve  optimum  performance,  it  is  necessary  to  understand  the  physics 
of  the  flow  structure  of  the  junction  vortex. 

A  number  of  researchers  have  studied  this  flow  using  numerical  and  experi¬ 
mental  methods.  Numerical  simulation  of  the  formation  of  the  junction  flow  is 
not  difficult  to  achieve  (Rubel,  1978)  since  the  generation  of  this  secondary  flow 
is  an  inviscid  process.  However,  downstream  of  separation,  a  mathematical  model 
becomes  complicated  because  of  a  lack  of  the  knowledge  of  the  mechanism  which 
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causes  the  jimction  vortex  phenomenon.  Recently,  computational  methods  for  the 
laminar  horseshoe  vortex  flow  were  developed  (Kaul,  Kwak,  and  Wagner,  1985; 
Briley  and  McDonald,  1985).  Comparison  with  the  experimental  data  in  the  sep¬ 
arated  zone  was  made  using  the  oil  flow  visualization  studies  and  wall  pressure 
measurements.  More  recently,  the  laminar  flow  fields  around  the  juncture  formed 
by  a  cylinder/flat  plate  were  simulated  using  the  full  three-dimensional  unsteady 
Navier-Stokes  equations  (Visbal,  1991).  Most  of  experimental  work  includes  flow 
visualization,  hot-wire  anemometer,  and  LDV  measurements. 

Flow  visualization  techniques  have  contributed  to  improving  the  physical  un¬ 
derstanding  of  complicated  flow  phenomena  (Kline,  1978),  and  have  been  an  es¬ 
sential  element  in  current  experimental  fluid  mechanics  research.  Especially  in 
turbulence  research,  flow  visualization  has  been  used  as  an  experimental  tool  to  in¬ 
vestigate  flow  structures  by  observing  complex  physical  phenomena  and  identifying 
them  as  combinations  of  simple  ones. 

In  most  of  the  studies  made  on  a  junction  flow,  a  surface  oil  flow  technique 
has  been  used  to  visualize  the  meein  flow  characteristics  on  the  surface  around  an 
obstacle  (Belik,  1973;  Sedney  and  Kitchens,  1975;  Baker,  1980;  Rood,  1984;  Moore 
and  Forlini,  1984;  Abid  and  Schmitt,  1986;  Eckerle  and  Langston,  1986;  Devenport 
and  Simpson,  1986,  1987,  1988,  1990;  Oilmen,  1990;  and  many  others).  It  has 
been  observed  that  there  is  a  system  of  vortices  (Sedney  and  Kitchens,  1975;  Baker, 
1979, 1980,  etc.).  The  flow  field  with  multiple  vortex  systems  was  constructed  using 
topology  and  some  flow  field  information  was  obtained  using  a  smoke  wire  visual¬ 
ization  and  local  fiow  field  measurements  (see  Figure  1).  However,  this  technique  is 
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confined  to  the  visualization  of  the  surface  region,  and  cannot  provide  the  charac- 
teristics  of  the  instantaneous  flow  field.  In  order  to  understand  the  flow  structure  in 
this  region,  a  method  is  required  that  can  visualize  three-dimensional  inotantaneous 
motion  in  the  flow. 

Many  researchers  have  studied  the  horseshoe  vortex  system  formed  around  a 
junction  using  various  measurement  techniques:  hot-wire  anemometers  (Shabaka 
and  Bradshaw,  1981;  McMahon  ti  al.,  1982;  Mehta,  1984;  Kubendran,  McMahon, 
and  Hubbartt,  1984;  Eckerle  and  Langston,  1986;  Devenport  and  Simpson,  1986, 
1990:  Greco,  1990),  pressure  measurements  (Baker,  1979,  1980;  Moore  and  For- 
lihi,  1984;  Eckerle  and  Langston,  1986;  Devenport  and  Simpson,  1986),  and  LDV 
measurements  (Tree,  1986;  Abid  and  Schmitt,  1986;  Devenport  and  Simpson,  1987, 
1990).  These  methods  were  .point  mcMurements  requiring  repetitive  measurements 
to  obtain  the-mean  flow  characteristics.  Point  by  point  measurements  are  necessary 
to  map  the  entire  mean  flow  field. 

Techniques  to  make  measurements  of  the  instantaneous  velocity  field  have  been 
developed  .with  the  aid  of  recent  advances  in  electronic  and  optical  technology.  Using 
a  photographic  method  with  image  processing,  the  velocity  field  can  be  mapped  by 
tracing  the  particles  present  in  the  photographic  images.  PIDV  (Particle  Image 
Displacement  Velocimetry)  is  the  method  that  has  been  used  widely  to  obtain  the 
velocity  field. 

In  this  research,  various  experimental  techniques  were  used  to  make  the  globed 
and  detailed  investigations  of  the  flow  structure  formed  around  the  junction.  The 
junction  flow  was  produced  using  a  flat  plate  and  a  wing  mounted  on  it.  The  wing 
was  composed  of  a  3  :  2  ellipticed  nose  an  NACA  0020  airfoil  tail  connected  at 
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its  maximum  thickness.  The  tests  were  conducted  in  a  water  tunnel  at  low  speeds. 
First,  hot-film  measurements  were  made  to  calibrate  the  flow  field  near  the  entrance 
of  the  test  section  of  the  water  tunnel.  LDV  measurements  were  performed  to  obtain 
the  averaged  velocity  in  the  plane  of  symmetry  upstream  of  the  wing.  The  results 
were  compared  with  the  wind  tunnel  data  (Devenport  and  Simpson,  1987)  measured 
on  a  half-size  wing  of  the  same  geometry  at  different  speed.  The  LDV  signals 
obtained  in  the  strong  vortex  region  were  analyzed  to  understand  the  low  frequency 
fluctuations  of  the  flow.  Emphasis  was  placed  on  developing  the  flow  visualization 
using  a  hydrogen  bubble  technique  and  the  application  of  the  PIDV  technique  using 
image  processing  techniques  to  extract  the  three-dimensional  instantaneous  velocity 
field  upstream  of  the  wing,  An  attempt  to  describe  the  flow  field  in  the  junction 
region  were  made  using  the  observations  obtained  from  the  flow  visualization. 

The  objectives  of  the  current  research  are  as  follows: 

(a)  to  visualize  the  flow  field  using  a  hydrogen  bubble  wire  technique  to  confirm 
the  previous  observation  of  the  multiple  vortex  system  (Sedney  and  Kitchens, 
1975;  Baker,  1979,  1980,  etc.), 

(b)  to  determine  the  connection  between  the  multiple  vortices,  backflow,  and  vortex 
stretching  using  the  flow  visualization  and  the  local  LDV  measurements, 

(c)  to  check  the  dependency  of  the  junction  vortex  behavior  on  the  Reynolds  num¬ 
ber, 

(d)  to  analyze  the  low  frequency  fluctuations  in  LDV  signals  to  obtain  the  physi¬ 
cal  understanding  of  the  bimodal  structure  found  by  Devenport  and  Simpson 
(1987,  1988a,  1988b,  1990), 

(e)  and  to  develop  the  PIDV  technique  to  obtain  the  3-D  instantaneous  velocity 
vector  field. 
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2.  BACKGROUND 


2.1  Investigations  of  a  Junction  Flow 

Previous  investigations  of  a  junction  flow  can  be  divided  into  two  main  cate¬ 
gories.  One  is  the  study  of  a  cylinder-end  wall  junction  flow  and  the  other  is  the 
study  of  the  idealised  wing-body  jimction  flow.  Both  experimental  models  use  a 
simplifled  geometry  to  reduce  the  complexity  of  the  flow;  a  flat  plate  for  the  devel¬ 
opment  of  the  oncoming  boundary  layer,  and  an  obstacle  (cylinder  or  wing  model) 
mounted  normal  to  it  for  deflecting  the  boundary  layer.  The  study  can  be  again 
categorized  according  to  the  flow  state  of  a  laminar  boundary  layer  or  a  turbulent 
boundary  layer. 


B^lik  (1973)  conducted  experimental  studies  of  Ihe  secondary  flow  about  a 
circular  cylinder  mounted  normal  to  a  flat  plate  to  simulate  a  flow  formed  in  blade 
cascades  bounded  by  two  parallel  flat  walls.  He  used  surface  oil  visualization  and 
static  pressure  measurements  on  the  flat  wall  in  the  separation  region  to  determine 
the  position  and  size  of  the  vortex  core.  He  obtained  r/il  w  1.43,  where  R  is  the 
radius  of  the  circular  cylinder  and  r  is  the  radial  distance  from  the  center  of  the 
obstacle  to  the  core  of  the  vortex.  This  was  in  good  agreement  with  the  estimates 
of  location  from  his  visualized  flow  pattern.  His  results  show  that  the  separation 
region  can  be  characterized  by  a  single  similarity  number  (the  vortex  Fourier  number 
{Fov  =  0-5p^)  or  the  Reynolds  number  based  on  cylinder  diameter  D). 
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Sedney  and  Kitchens  (1975)  performed  investigations  of  the  structure  of  three- 
dimensional  ^parated  flows  in  obstacle,  boundary  layer  interactions.  Their  study 
was  conducted  in  a  turbulent,  supersonic  boundary  layer  flows  (M  =  1.5,  2.5,  3.5, 
4.0,  and  4.5)  over  protuberances.  They  used  small  obstacles,  which  had  k  <  6  where 
k  is  the  obstacle  height  and  S  is  the  99  %  boundary  layer  thickness.  They  found 
that  the  cylinder  height  affects  the  strength  and  characteristics  of  the  vortex  when 
k<S. 

Baker  (1979)  investigated  the  horseshoe  vortex  formed  around  the  base  of  a 
cylinder  of  diameter  D  by  a  separating  laminar  boundary  layer.  He  made  a  dimen¬ 
sional  analysis  on  (the  distance  from  the  obstacle  center  to  the  center  of  vortex 
core),  X,  (the  distance  from  the  obstacle  center  to  the  separation  line  upstream 
of  the  obstacle  in  the  plane  of  symmetry),  and  /  (the  frequency  of  oscillations  of 
the  horseshoe  vortex  system),  and  the  corresponding  dimensionless  group  can  be 
expressed  as  follows: 


where 


£l  l£.  -  £.\ 

D'  U  u  ' 

is  the  Strouhal  number. 


(2.1) 


Using  a  smoke-wire  visualization  technique,  he  observed  three  different  types  of  vor¬ 
tex  systems,  which  were  ;  i)  steady  horseshoe  vortex  systems  with  2,  4  or  6  vortices, 
the  number  of  vortices  increasing  as  UDfv  increases,  ii)  horseshoe  vortex  systems 
which  exhibit  a  regular  oscillatory  motion,  and  u'»)  horseshoe  vortex  systems  which 
exhibit  an  irregular  unsteady  behavior.  Spectral  analyses  of  the  hot-wire  signals 
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were  implemented  to  investigate  the  oscillating  behavior  of  the  horseshoe  vortex 
system.  He  made  explanations  that  at  certain  values  of  UD/v  and  U/S*  the  horse¬ 
shoe  vortex  itself  becomes  unstable  for  some  reason,  and  begins  to  oscillate  in  one 
of  two  ‘natural’  modes,  the  higher  frequency  mode  becoming  more  dominant  as  the 
Reynolds  number  increases.  He  found  that  Stj}  =  0.26  for  low  frequency  oscillation, 
and  Sto  —  0.40  for  high  frequency  oscillation  (5<d  increased  up  to  0.60  at  higher 
UDju) 

For  the  turbulent  horseshoe  vortex,  Baker  (1980)  made  experimental  investi¬ 
gations  of  the  model  of  the  same  geometry.  He  observed,  by  the  surface  oil-flow  vi¬ 
sualization,  that  only  the  vortex  system  consisting  of  four  vortices  (see  Figure  1(b)) 
exists  for  the  test  conditions  used  (  4000  <  VD/v  <  90000,  4  <  DjS*  <  30  ). 

Eckerle  and  Langston  (1986)  made  measurements  of  a  turbulent  horseshoe  vor¬ 
tex  formed  around  a  cylinder.  They  made  detailed  surface  flow  visualization  and 
static  pressure  measurements,  and  extensive  mean  velocity  and  pressure  measure¬ 
ments  in  and  around  the  flow  fleld  of  the  vortex  system.  They  observed  only  a 
single  horseshoe  vortex  in  their  visualization. 

The  study  by  Greco  (1990)  showed  the  flow  behavior  near  a  cylinder-end  wall 
junction  using  hydrogen  bubble  flow  visualization  and  hot-film  anemometry  in  a 
low-speed,  free-surfacs  water  channel.  Five  distinct  laminar  flow  regimes  were 
investigated,  and  he  documented  that,  in  a  periodic  flow  regime,  the  periodicity 
of  the  generation  of  the  junction  vortex  was  found  to  be  a  function  of  free-stream 
velocity,  cylinder  diameter,  fluid  viscosity,  and  streamwise  location.  He  observed 
that  above  a  critical  Reynolds  number,  the  vortex  system  breaks  down  to  become 
a  large-scale  turbulent  horseshoe  vortex. 


BACKGROUND 


8 


Idealized  wing-bodv  iufiction 

Schwind  (1962)  made  a  study  on  an  incompressible  laminar  boimdary  layer 
separation  at  the  junction  of  a  wall  and  a  60“  included-angle  wedge.  He  used 
smoke-wire  visualization  and  observed  five  flow  regimes  according  to  the  flow  char¬ 
acteristics;  three  steady  flow  regimes  at  low  speeds,  and  two  periodic  flow  regimes 
at  higher  flow  velocities. 

Shabaka  and  Bradshaw  (1981)  made  extensive  cross-wire  measurements  of  tur¬ 
bulent  quantities  in  an  idealized  wing-body  junction.  The  wing  consisted  of  a  semi¬ 
elliptical  leading  edge  followed  by  a  slab  of  constant  thickness.  Their  measurements 
were  confined  to  the  region  of  nearly  constant  pressure  downstream  of  the  leading 
edge.  The  results  showed  that  the  vortex  decayed  only  slowly  with  distance  down¬ 
stream,  and  its  diameter  remained  roughly  equal  to  the  thickness  of  the  boundary 
layer  on  the  body.  They  found  that  the  attenuation  of  the  skew-induced  secondary 
flow  by  the  Reynolds  stresses  was  very  slow. 

McMahon,  Hubbartt,  and  Kubendran  (1982)  measured  three  mean  velocity 
components  and  six  turbulent  stresses  in  a  juncture  flow  formed  by  a  constant 
thickness  body  having  an  elliptical  leading  edge  (which  was  mounted  perpendicu¬ 
lar  to  a  large  flat  plate).  Two  single  hot  wire  sensors  were  used  at  two  different 
streamwise  stations.  They  found  that  the  strength  of  the  secondary  flow  vortex  in 
the  juncture  increases  as  the  leading  edge  of  a  body  of  constant  thickness  is  made 
more  blunt,  and  that  this  stronger  vortex  has  an  effective  core  which  is  closer  to 
the  surface  of  the  body  than  that  of  a  weaker  vortex  caused  by  a  leading  edge  of  a 
smellier  flneness  ratio. 
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Mehta  (1984)  studied  the  effects  of  the  wing  nose  shape  on  the  flow  in  a  wing- 
body  junction.  He  used  the  same  model  as  Shabaka  and  Bradshaw  (1981)  except 
that  he  had  three  nose  shapes.  Using  miniature  cross-wire  probes,  the  secondary 
flow  was  measured  while  the  streamwise  velocities  were  measured  using  standard 
pitot  tubes.  It  was  found  that  the  nature  of  the  horseshoe  vortex  is  strongly  depen¬ 
dent  on  the  wing  nose  shape.  The  size  and  the  strength  of  the  vortex  increase  with 
nose  bluntneas. 

Moore  and  Forlini  (1984)  investigated  the  horseshoe  vortex  formed  around  a 
Rankine  half-body.  The  body  was  located  in  a  duct  bounded  by  two  plane  endwalls 
and  two  side  walls  shaped  like  potential-flow  streamlines.  They  made  surface  oil- 
flow  visualization  and  static  pressure  measurements  on  the  endwall  and  on  the 
’  surface  of  the  body  on  the  planes  downstream  of  the  leading  edge.  A  five-hole 
pressure  probe  was  traversed  to  measure  the  secondary  flow  in  their  measurements. 

Abid  and  Schmitt  (1986)  used  a  directionedly  sensitive,  three-component  laser 
Doppler  velocimeter  to  make  the  mean  velocities  and  turbulent  Reynolds  measure¬ 
ments  with  oil  flow  visualization  and  static  pressure  measurements  at  the  wall.  The 
horseshoe  vortex  system  was  generated  by  the  obstacle  mounted  perpendicular  to 
a  flat  plate.  The  leading  edge  of  the  obstacle  was  a  semi-circular  cylinder  and  its 
downstream  part  had  a  streamlined  shape  intended  to  avoid  additional  flow  separa¬ 
tion  and  to  keep  the  flow  as  steady  as  possible.  They  presented  LDV  measurement 
data  in  the  separated  region. 

Kubendran,  McMahon,  and  Hubbartt  (1986)  made  measurements  of  mean  ve¬ 
locity  components  and  turbulent  stresses  in  the  entire  flow  fleld  surrounding  the 
leading  edge  of  a  wing/fuselage-type  juncture.  A  constant  thickness  body  with  a 
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1.5  :  1  leading  edge  was  mounted  normal  to  a  flat  plate.  They  found  that  both  the 
strength  and  the  location  of  the  horseshoe  separation  vortex  can  be  controlled  by 
modifying  the  leading  edge  shape  of  the  wing. 

Dickinson  (1986a,  1986b)  presented  LDV  measurements  for  the  geometrically 
similar  model  as  the  one  used  in  this  research.  He  reported  that  the  larger  cross- 
flow  velocities  resulted  with  increased  bluntness  of  the  nose  shape.  He  used  oil 
film  and  oil  dot  flow  visualization  techniques  (1988)  to  find  surface  separation  and 
attachment  lines  and  streamline  patterns. 

Devenport  and  Simpson  (1986,  1987,  1988a,  1988b,  1990)  have  contributed 
to  the  understanding  of  the  time-dependent  behavior  of  the  turbulent  wing-body 
junction  flow  though  their  experiments.  Using  a  single  hot-wire  probe  and  a  three 
component  laser  doppler  anemometer,  they  made  measurements  of  the  mean  ve¬ 
locity  and  turbulent  stresses  in  the  region  surrounding  the  wing.  Instantaneous 
spanwise  distributions  of  the  skin  friction  were  obtained  using  an  array  of  surface 
hot-wire  sensors.  Surface  pressure  measurements  and  oil  flow  visualization  were 
conducted  on  the  flat  wall  surrounding  the  wing.  They  found  bimodal  velocity 
histograms  in  the  junction  region,  related  to  low  frequency  fluctuations  in  the  flow. 

Devenport  et  al.  (1989)  made  mean  pressure  and  hot-wire  measurements  to 
show  the  effects  of  a  simple  fillet  on  the  flow  of  a  turbulent  boundary  layer  past  an 
idealized  wing-body  junction. 
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2.2  Hydrogen  Bubble  Technique 


Hydrogen  bubbles  are  used  as  tracer  particles  to  visualize  the  flow.  The  hydro¬ 
gen  bubble  method  is  based  on  the  electrolysis  of  water.  Although  the  ions  depend 
upon  the  electrolyte  involved,  in  this  case,  the  reaction  is  the  decomposition  of 
water  {2H2O  -*  2H2  +  O2).  When  two  electrodes  with  a  potential  difference  are 
immersed  in  water,  hydrogen  is  generated  at  the  cathode  and  oxygen  at  the  anode. 
If  the  cathode  is  a  thin  wire,  very  fine  hydrogen  bubbles  are  formed  —  the  bubbles 
released  from  the  wire  are  one-half  the  wire  diameter.  A  wire  may  be  used  as  either 
an  anode  or  a  cathode,  but  since  the  hydrogen  bubbles  are  produced  at  twice  the 
rate,  they  are  used  as  flow  tracers  in  this  technique. 

The  hydrogen  bubble  technique  was  first  described  in  1953  by  Kolins.  Using 
a  copper  wire  across  a  Lucite  tube,  he  visualized  the  parabolic  distribution  of  the 
velocity  in  a  laminar  pipe  flow.  A  fundamental  improvement  ir  this  method  was 
discovered  by  Geller  (1955)  when  he  investigated  the  boundary  layer  on  a  water 
tunnel  wall.  He  began  to  use  a  thin  platinum  wire  of  diameter  of  0.025mm  as  a 
cathode  in  his  study.  This  technique  was  further  developed  by  Clutter  and  Smith 
(1961). 

In  the  1960s,  this  technique  began  to  prevail  in  investigating  the  turbulent 
flows  because  of  its  ease  and  flexibility  in  application.  Schraub  et  al.  (1965)  applied 
the  hydrogen  bubble  technique  to  the  study  of  time-dependent  low  speed  water 
flows.  They  presented  detailed  descriptions  of  various  applications  of  this  technique. 
In  particular,  a  time-line  and  a  combined  time-streak  marker  method  provided 
quantitative  information  on  the  time-dependent  velocity  and  trajectories  of  the 
fluid  elements. 
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Asaouma  and  Takeda  (1965)  and  Mattingly  (1966)  extended  the  use  of  this 
technique  to  higher  velocities  (up  to  7m/s  and  5m/s,  respectively)  using  a  higher 
cathode  potential  and  increased  electrolyte  additive  concentration. 

Criminale  and  Nowell  (1965)  showed  that  this  technique  could  be  extended  to 
the  measurements  in  liquid-water  mixtures.  They  succeeded  in  using  this  method 
in  a  glycerine-water  mixture  up  to  40  %  (by  weight)  of  glycerine.  Roos  and  Will- 
marth  (1969)  reported  some  results  obtained  through  use  of  the  hydrogen  bubble 
technique  in  high  concentration  glycerine  mixtures  in  their  study  of  low  Reynolds 
number  unsteady  wake  flows.  In  particular,  they  provided  the  justification  for  using 
hydrogen  bubbles  as  tracer  particles  by  showing  that  the  relative  motion  between 
the  bubble  and  the  surrounding  fiuid  is  negligible. 

Thompson  reviewed  this  technique  in  his  paper  (1972):  His  discussion  included 
the  number  of  bubbles  produced  at  the  cathode,  the  variation  of  cathode  current 
with  electrode  potential,  variation  of  water  path  resistance  with  electrolyte  concen¬ 
tration,  effect  of  anode  size  and  electrode  spacing  on  water  path  resistance,  and  the 
effect  water  temperature  on  water  path  resistance. 

Kim  ei  al.  (1968,  1971)  and  Grass  (1971)  used  this  technique  to  study  the 
turbulent  boundary  layer  near  smooth  and  rough  walls,  respectively.  McAllister 
and  Carr  (1978)  investigated  the  dynamic  stall  over  an  oscillating  airfoil  using  this 
technique.  In  the  latter  study,  a  wire  cathode  was  stretched  across  the  test  section 
and  oriented  normal  to  the  direction  of  the  flow. 

Smith  and  Paxspn  (1983)  at  Lehigh  University  developed  the  computer-augmented 
hydrogen  bubble-wire  flow  visualization  system.  Using  digitally  interfaced  dual¬ 
view  video  sequences  of  a  bubble  time-line,  the  three  dimensional  motion  of  the 
time-line  was  recreated,  and  the  Lagrangian  velocity  was  obtained. 
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Appli<»tions  of  this  technique  has  greatly  increased  since  its  introduction  due  to 
increasing  need  to  understand  global  flow  patterns  and  improvements  in  electronic 
wd  optical  instrumentation.  In  a  turbulent  flow,  the  hydrogen  bubble  method  has 
played  an  important  role  as  an  experimental  tool  for  identifying  new  flow  structures. 

2.3  Image  Processing  Technique  Applications 

Flow  visualization  resulting  from  the  interaction  between  light  and  matter  has 
been  useful  to  design  experiments  and  for  proper  interpretation  of  the  conventional 
measurements.  Recently,  rapid  advances  in  scientific  technology  have  enabled  flow 
visualization  to  be  used  as  a  direct  measurement  tool  for  the  flow  characteristics 
through  the  application  of  digital  image  processing  techniques. 

Our  motivation  for  using  digital  image  processing  was  the  current  emphasis  on 
detecting  and  classifying  large-scale  structures  in  flows.  Employing  digital  picture 
processing  techniques,  the  information  of  the  entire  flow  field  frozen  in  the  flow 
visualization  pictures  can  be  measured,  while  conventional  measurement  techniques 
such  as  hot-wire  anemometry  and  LDV  measurements  only  allow  single  point  mea¬ 
surements. 

One  application  of  digital  image  processing  to  flow  visualization  is  particle  im¬ 
age  velocimetry  (PIV)  or  particle  image  displacement  velocimetry  (PIDV).  Accord¬ 
ing  to  the  classification  of  digital  image  processing  and  flow  visualization  techniques 
by  Hesselink  (1988),  this  belongs  to  the  group  that  obtains  the  images  by  scattering. 
Two  other  distinct  approaches  that  have  been  widely  applied  using  scattering  are 
laser  speckle  velocimetry  (LSV),  and  laser  induced  fluorescence  (LIF). 
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This  group  uses  light  sheet  generated  by  a  laser  beam  to  scatter  light  from  the 
seeding  particles  in  the  flow.  Since  PIDV  uses  a  low  density  of  seeding  particles, 
the  measurements  are  possible  only  at  those  locations  where  particles  happen  to  be 
during  exposure.  Difficulties  that  arise  in  the  application  of  PIDV  are  that  parti¬ 
cles  are  indistinguishable  from  each  other  unless  special  care  is  taken  and  particles 
moving  out  of  the  illumination  thickness  of  the  light  sheet  may  not  be  distinguished 
from  the  particles  staying  inside.  The  major  disadvantage  of  PIDV  is  that  this 
technique  can  be  applied  only  to  the  measurements  in  one  plane  so  that  it  is  not 
easy  to  use  the  PIDV  in  regions  where  one  is  interested  in  three-dimensional  flow. 

A  number  of  researchers  have  used  the  flow  visualization  pictures  to  obtain  the 
quantitative  information  of  the  entire  flow  field.  Grass  (1971)  used  hydrogen  bubble 
flow  tracers  and  medium  high-speed  motion  photography  to  obtain  an  instantaneous 
visual  and  quantitative  description  of  the  flow  field.  Using  instantaneous  longitudi¬ 
nal  and  vertical  velocity  profiles,  he  computed  the  distributions  of  the  instantaneous 
Reynolds  stress  contribution.  Imaichi  and  Ohmi  (1983)  developed  a  system  for  eval¬ 
uating  some  of  the  physical  variables  of  fluid  flows  through  flow  visualization  and 
image  processing  techniques.  Their  measurements  were  within  reasonable  accuracy 
for  two-dimensional  vortical  flow  examples. 

Utami  and  Ueno  (1984)  developed  a  tracer  method  to  visualize  the  three- 
dimensional  flow  structure  of  a  turbulent  open  channel  flow.  They  designed  a 
movable  lighting  system  in  which  horizontal  cross-sections  of  the  flow  could  be 
illuminated  instantaneously  in  a  sequence  by  shifting  the  lighting  plane  using  a 
stepping  motor.  A  Lagrangian  description  of  the  flow  field  was  used  in  the  presen¬ 
tation  of  their  results. 
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Lu  and  Smith  (1985)  employed  automated  image  processing  of  hydrogen  bubble 
flow  visualization  pictures  to  establish  local,  instanteineous  velocity  profile  informa¬ 
tion.  They  used  the  time-of-flight  technique  to  construct  the  local  velocity  profiles 
from  the  hydrogen  bubble  images,  and  attempted  to  correlate  prototype  turbulent 
burst  detection  techniques  with  corresponding  visualization  data. 

Adamczyk  and  Rimm  (1988a)  used  two-dimensional  particle  tracking  velocime- 
try  (PTV)  to  determine  the  fluid  velocity  in  a  transparent  test  section  simulat¬ 
ing  an  engine  throttle  body  assembly.  They  further,  developed  automated  three- 
dimensional  PTV  using  two  cameras  in  their  following  study  (Adamczyk  and  Rimai, 
1988b). 

Racca  and  Dewey  (1988)  developed  a  technique  to  automatically  analyze  the 
three-dimensional  motion  of  tracers  in  a  fluid  flow.  They  designed  the  split  field 
mirror  system  in  ordei  to  obtain  simultaneous  orthogonal  views  in  a  single  high 
speed  cine  camera.  Eulerian  information  weis  derived  from  the  pooled  velocity  data 
points  by  interpolation  on  a  regular  spatial  grid. 

Ensemble-averaging  and  correlation  techniques  were  developed  by  Kerstens 
and  Rockwell  (1988)  when  they  studied  macroscale  features  of  unsteady  flows.  They 
used  hydrogen  bubble  time-line  images  for  tracking  and  averaging. 

Liu,  Landreth,  Adrian,  and  Hanratty  (1991),  improved  the  PIDV  technique 
by  developing  high  resolution  particle  image  velocimetry  to  obtain  the  turbulent 
velocity  field  for  fully  developed  flow.  They  made  comparison  of  their  PIDV  results 
with  those  of  direct  numerical  simulation  and  LDV  measurements,  and  showed  good 
agreement. 
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3.  EXPERIMENTAL  SETUP  AND  APPARATUS 


3.1  Water  Tunnel 

The  water  timnel  of  the  Aerospace  and  Ocean  Engineering  Department  at 
Virginia  Polytechnic  institute  and  State  University  is  a  horizontal  closed-loop  type 
with  test  section  size  of  0.61m  x  0.61m  x  2.44m.  The  water  timnel  is  made  up  of 
circular  cross  sections  of  0.46m  O.D,  PVC  pipe  except  the  rectangular  honeycomb- 
screen  section  and  test  section.  The  flow  is  driven  by  an  Ingersoll-Rand  0.46m  I.D. 
axial  pump  (Model  No.  B0104DLF2yD)  connected  through  a  timing  belt  to  a  3 
phase,  10  hp  induction  motor,  which  is  controlled  by  a  Toshiba  TOSVERT-130 
H/Hl  low  acoustic-noise  transistorized  inverter. 

Three  honeycomb  and  eight  screen  sections  are  located  upstream  of  the  test 
section  to  isolate  the  test  section  from  the  pump  and  improve  the  uniformity  of  the 
flow  (see  Hsu  (>1989)  for  detailed  description).  The  test  section  is  made  of  1.27cm 
thick  glass  to  provide  a  visual  access  to  the  working  area.  Figure  2  shows  the  general 
layout  of  the  water  tunnel. 

The  flow  in  the  test  section  was  surveyed  for  two  speed  cases  using  hot  film 
measurements  at  three  stations  with  the  wing-body  junction  model  installed.  The 
flow  velocities  surveyed  were  21cm /see  and  49cm/sec,  respectively.  The  turbulence 
level  of  the  low  speed  case  at  the  entrance  of  the  test  section  ranged  from  0.4% 
to  1.6%  in  the  plane  of  symmetry.  For  the  high  speed  case,  the  turbulence  level 
of  0.8%  -  1.3%  was  obtained  at  the  same  location.  It  was  observed  that  as  the 
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wing  is  approached  the  turbulence  level  becomes  higher  in  the  lower  center  which  is 
closest  to  the  wall  among  the  measuring  stations.  The  explanation  of  the  increasing 
turbulence  level  as  the  flow  approached  the  wing  may  be  made  as  due  to  several 
possible  reasons,  which  may  be  the  effect  of  the  axial  pump  not  reduced  even  through 
the  honeycomb  section  or  impurities  in  the  screens.  However,  it  is  believed  that  this 
does  not  influence  the  experimental  results  and  conclusions  obtained  in  this  study. 
The  measuring  locations  and  the  turbulence  level  at  each  location  are  as  shown  in 
Figure  3. 

3.2  Model 

The  wing-body  junction  model  consists  of  a  wing,  a  false  floor,  and  two  side 
walls,  all  made  of  plexiglas  for  the  purpose  of  flow  visualization.  The  assembled 
wing-body  junction  model  mounted  in  the  test  section  is  as  in  Figure  4. 

3,2,1  Wing 

The  wing  hris  a  61cm  chord,  a  58.4cm  height,  and  a  14.35cm  thickness.  The 
cross  section  of  the  wing  is  composed  of  a  3  :  2  elliptical  nose  aind  an  NACA  0020 
airfoil  section,  which  are  connected  at  the  maximum  thickness  location.  Plexiglas 
was  used  to  build  the  wing  model  to  provide  a  view  looking  upstream  at  the  flow 
in  the  nose  region  of  wing-body  junction.  The  view  was  accessed  from  the  bottom 
through  the  wing,  where  a  45°  inclined  plane  mirror  weis  installed  to  face  upstreeim. 
The  wing  was  reinforced  using  five  inner  aluminum  ribs  coated  with  marine  epoxy. 

A  thermal  forming  technique  weis  employed  to  obtain  the  leading  edge  shape  of 
the  plexiglas  wing.  A  plywood  form  was  so  built  that  the  thickness  of  the  plexiglas 
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(3.175mm)  was  taken  into  account  for  the  outermost  contour  to  have  the  desired 
shape.  The  form  and  a  square  of  3.175mm  thick  plexiglas  placed  over  it  were  heated 
to  310®  F  in  an  oven  until  the  softened  plexiglas  reached  the  proper  shape.  After 
the  plexiglas  obtained  approximately  correct  contour,  minor  manual  adjustment 
was  required  to  trim  the  desired  shape.  Natvural  cooling  in  the  oven  was  necessary 
until  a  room  temperature  was  reached.  When  cooled,  the  plexiglas  was  deformed 
along  the  outer  edge.  Therefore,  the  pieces  were  made  oversized  so  that  they  could 
be  trimmed  after  cooling. 

5.2.2  False  Floor 

A  false  floor  consisted  of  a  horizontal  flat  plate  designed  to  mount  the  wing  and 
join  two  side  walls.  It  was  constructed  to  maintain  visual  access  from  the  bottom  of 
the  test  section.  This  floor  started  after  the  contraction  section.  The  total  length  of 
the  false  floor  from  tho  entrance  to  the  leading  edge  of  the  wing  was  approximately 
1.19m. 

5.2.3  Side  Walls 

To  compensate  for  the  blockage  effect  due  to  the  large  size  of  the  wing  model, 
two  flexible  side  walls  were  mounted  vertically  on  the  false  floor.  The  outer  contour 
of  the  side  walls  was  shaped  to  match  the  streamlines  calculated  by  an  inviscid 
source  panel  method  with  a  boundary  layer  displacement  thickness  correction.  To 
determine  the  final  contour  of  the  side  walls,  the  displacement  thickness  was  cal¬ 
culated  using  a  turbulent  flat  plate  boundary  layer  originated  from  the  end  of  the 
contraction  section.  This  values  of  displacement  thickness  were  added  to  the  Y 
locations  of  streamline  obtained  by  the  potential  flow  calculation.  The  location  of 
the  side  walls  with  and  without  a  displacement  correction  is  compared  in  Figure  5. 
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S.i.4  Contraction  section 

A  1.6:1  contraction  section  was  built  to  prevent  the  flow  separation  which 
otherwise  would  occur  at  the  leading  edge  of  the  test  section  floor  and  two  side 
walls.  The  section  was  located  immediately  downstream  of  the  last  screen  and 
connected  at  the  leading  edge  of  the  false  floor. 
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4.  EXPERIMENTAL  TECHNIQUES 


4.1  Hydrogen  Bubble  Technique 

4.1.1  General  Description 

A  0.025mm  tungsten  wire  was  employed  as  the  cathode  immersed  in  the  test 
section  of  the  water  tunnel.  The  anode  was  a  piece  of  aluminum  foil  and  floated 
aft  of  the  wing  model  to  avoid  disturbance  of  the  flow.  We  used  a  low  speed 
{U  =  9.07cm/s  producing  Reg  =  330)  to  prevent  cavitation-induced  air  bubbles 
which  interfere  with  the  small  hydrogen  bubbles.  Hydrogen  bubbles  produced  at  the 
wire  were  swept  off  by  the  flow  and  followed  downstream.  The  predominant  force 
on  the  hydrogen  bubbles  was  the  drag  due  to  the  local  water  motion.  Generally,  the 
buoyancy  force  on  the  bubbles  is  greatly  reduced  by  decreasing  the  wire  diameter, 
and  was  almost  negligible  at  the  wire  diameter  of  0.025mm  (Schraub  et  al.  1965). 
In  this  situation,  the  hydrogen  bubbles  can  be  considered  to  follow  the  local  flow 
motion.  It  is  noted  that  the  velocity  defect  due  to  the  wire  should  be  taken  into 
account  before  the  bubbles  reach  the  local  flow  speed.  According  to  Schraub  et  al. 
(1965),  The  bubble  velocity  reaches  the  free-stream  velocity  in  less  than  70  wire 
diameter,  which  corresponds  to  1.75mm  in  this  study.  A  small  amount  of  additive 
(  0.15  gr  of  Na2S04  per  a  liter  of  water  as  suggested  by  Schraub  et  al.  (1965))  was 
used  to  enhance  the  contrast  with  the  increeised  production  of  hydrogen  bubbles. 

In  this  study,  we  developed  three  techniques  as  we  applied  the  hydrogen  bubble 
flow  visualization.  First  of  all,  a  time-line  technique  was  used  to  visualize  the  flow 
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structure  near  the  nose  region  of  the  wing  model.  The  time-line  patterns  were 
obtained  by  producing  the  potential  difference  between  the  wire  and  the  anode 
in  square-wave  pulses.  The  frequency  of  the  time-lines  was  adjustable  through 
the  function  generator  connected  to  the  hydrogen  bubble  generator  circuit.  The 
hydrogen  bubble  probe  for  the  time-line  patterns  had  a  vertical  wire  of  7.62cm 
height  and  a  support  to  hold  the  wire  in  a  vertical  position  in  the  symmetry  plane. 
The  details  of  the  probe  are  described  in  the  next  section. 

A  continuous  sheet  technique  is  analogous  to  the  smoke-wire  technique  used 
in  the  wind  tunnel  flow  visualization.  The  continuous  sheet  of  hydrogen  bubbles  is 
generated  by  supplying  the  circuit  with  a  constant  voltage  difference.  The  time-line 
and  a  continuous  sheet  technique  can  be  combined  to  generate  inverted  time-lines 
using  series  of  inverted  pulses.  This  technique  had  an  advantage  that  the  individual 
time-lines  could  be  seen  while  the  flow  field  was  still  observable  with  hydrogen 
bubbles  just  like  the  smoke-wire  technique. 

A  combined  time-streak  marker  technique  was  developed  to  get  quantitative 
data  from  flow  visualization  by  picture  processing.  Here,  the  bubble  generating  wire 
was  insulated  with  fine  enamel  at  regular  intervals  so  that  approximately  square 
areas  of  bubbles  were  produced  when  pulsed  by  the  generator  circuit  as  used  in  the 
time-line  technique.  A  brush  with  a  very  fine  tip  and  a  pattern  template  were  used 
to  mask  the  wire. 

A  multiple  wire  hydrogen  bubble  technique  was  developed  to  visualize  the 
three-dimensional  flow  structure  by  generating  bubble  patterns  by  the  sequenced 
pulsing  to  the  wires.  Since  the  probe  was  located  near  the  wing  and  flow  speed  at 
each  wire  was  different,  the  sequence  of  pulsing  to  each  wire  was  chosen  to  avoid 
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bubbles  overlapping  in  the  image  plane.  This  may  take  place  though  the  bubbles 
are  independent  in  three-dimensional  space.  The  details  of  the  probe  and  circuit 
design  Me  described  in  the  next  section.  The  sequenced  multiple  wire  hydrogen 
bubble  technique  used  the  combined  time-streak  marker  method.  Figure  6  shows 
an  example  of  the  sequence  of  hydrogen  bubble  generation  in  the  upstream  view 
with  the  corresponding  side  view.  Examples  of  each  application  are  presented  in 
Figure  7. 

4.I.B  Probe  and  Circuit  Design 

We  used  a  D-shaped  probe  with  width  of  22.86cm  and  height  of  17.78cm.  The 
probe  was  reinforced  with  cross  bars  and  it  was  connected  to  the  hydrogen  bubble 
generator  circuit  outside  the  test  section.  The  probe  designs  for  each  application 
are  illustrated  in  Figure  8. 

Each  cross-section  of  the  probe  had  a  streamlined  shape  to  minimize  flow 
interference  when  placed  in  the  flow.  A  Kevlar  holder  was  built  to  form  a  tight 
rectangular  loop  edong  the  edges  of  the  D-shaped  probe.  Four  glass  posts  of  length 
of  2.54cm  and  diameter  of  3.18mm  were  installed  at  two  lower  ends  and  8.89cm 
above  the  probe  leg.  These  were  used  to  make  a  stiff  electrically  isolated  holder. 

Two  probes  were  designed  for  each  application.  A  single  0.025mm  tungsten 
wire  was  moimted  vertically  in  the  center  of  the  probe  as  shown  in  Figure  8(b). 
Time-line  patterns,  continuous  sheet  and  inverted  time-lines  of  hydrogen  bubbles 
used  the  single  wire  probe. 

Five  0.025m77^  tungsten  wires  were  mounted  on  a  Kevlar  holder  with  a  1.27cm 
spacing  and  connected  to  the  top  of  the  probe.  All  but  8.89cm  of  wires  from  the 
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bottom  were  masked  with  enamel  in  order  not  to  generate  the  hydrogen  bubbles 
outside  the  region  of  interest.  This  probe  is  as  shown  in  Figure  8(c)  and  was  used 
for  the  sequenced  multiple  wire  technique. 

A  bubble  generator  circuit  was  designed  and  built  to  generate  the  pulses  at 
an  adjustable  time  constant  interval.  These  were  then  distributed  to  each  wire  in 
sequence.  The  interval  between  pulses  and  the  order  of  bubble  generating  wires 
were  adjustable.  Mercury  relays  were  used  to  reduce  the  bouncing  effect  of  the 
switch.  A  schematic  diagram  of  the  generator  circuit  is  as  shown  in  Figure  9. 

4>LS  Lighting 

Lighting  is  cruciai  in  the  hydrogen  bubble  technique.  The  lighting  setup  is 
described  in  Figure  10. 

For  still  photography,  we  used  two  flash  lights  located  at  opposite  side  of  the 
test  section  to  the  camera  to  pioduce  forward-scattered  light  from  the  hydrogen 
bubbles.  Two  synchronized  flashes  were  used  to  eliminate  the  shade  created  by  the 
probe.  The  location  of  the  flashes  was  restricted  by  the  external  configuration  of 
the  test  section,  such  as  the  threaded  rods  located  vertically  at  every  foot  outside 
the  test  section.  We  used  ASA  400  film  at  a  setting  of  l/125sec  and  an  f  number 
of  16  for  still  pictures. 

We  needed  more  light  for  the  high  speed  video  camera  than  for  still  photog¬ 
raphy.  For  a  side  view,  a  600W  quartz  halogen  lamp  was  located  outside  the  test 
section  at  an  angle  inclined  to  the  test  section  and  configured  to  yield  forward  scat¬ 
ter  on  the  hydrogen  bubbles.  A  second  lamp  was  needed  to  eliminate  shading  by 
the  probe.  For  a  bottom  view,  three  thousand-watt  focusable  halogen  lamps  were 
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used  for  the  illumination  of  the  hydrogen  bubbles  from  the  top  of  the  test  section 
(see  Figure  10).  Yet,  this  was  not  enough  to  get  high  quality  images  since  water 
absorbed  most  of  the  light  energy  on  the  way  to  the  bubble  region. 

Approximately  60®  inclined  to  the  line  normal  to  the  test  section  wall  produced 
an  effective  scattering  on  hydrogen  bubbles. 

4.2  Laser  Doppler  Velocimetry 

A  dual  beam,  forward-scattered,  frequency-shifted  LDV  (see  Figure  11)  was 
used  to  measure  the  axial  component  ve’coities  and  the  turbulent  quantities  in  the 
plane  of  symmetry  upstream  of  the  nose  of  the  wing-body  junction. 

4-2.1  General  Description  of  LDV 

Frequency  shifted  laser  velocimetry  was  first  introduced  by  Yeh  and  Cummins 
(1964).  They  added  an  artificial  frequency  difference  between  the  two  incident  laser 
beams  to  generate  a  set  of  moving  fringe  at  the  intersection  of  the  beams.  Since 
the  fringes  are  moving  in  one  direction,  even  a  stationary  particle  will  produce  a 
Doppler  frequency  at  the  shifted  frequency.  Consequently,  a  negative  velocity  of 
maximum  value  corresponding  to  less  than  the  shift  frequency  can  be  measured. 

When  a  particle  passes  through  the  intersection  volume  of  two  beams,  a  pedestal 
is  generated  along  with  the  Doppler  frequency.  A  pedestal  is  related  to  the  amount 
of  time  taken  for  a  particle  to  pass  through  the  probe  volume  and  is  removed  using 
a  high-pass  filter.  The  Doppler  frequency  produced  by  the  interference  pattern  has 
a  direct  linear  relationship  with  the  particle  velocity  as  follows. 
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U=M- 

2smK 
=  fodf. 

where  U  =  flow  velocity, 

fo  =  Doppler  frequency, 

A  =  wave  length  of  laser  beam, 

K  =  half  angle  between  two  incident  beams, 
and  df  =  fringe  spacing. 

With  a  frequency  shift,  Equation  (4.1)  is  modified  as  follows, 


where 


U  +  Uf,= 


ifn  +  /.)A 


2sinK 
=  {Id  +  ft)df 

ft  is  the  shifted  frequency. 


(4.1) 


(4.2) 


Using  Equation  (4.2),  the  flow  velocity  can  be  calculated  with  the  known  fringe 
spacing  and  shifted  frequency. 


4-2-2  LDV  Setup 

The  LDV  consisted  of  an  INNOVA  70  argon-ion  laser  operated  at  514.5  /im. 
As  shown  in  Figure  11,  the  beam  from  this  laser  was  reflected  using  severed  mirrors 
to  be  directed  to  the  test  section.  A  beam  splitter  was  used  to  form  two  incident 
beams.  After  passing  through  two  Bragg  cells,  two  approximately  40  Mhz  shifted 
first  order  beams  were  generated.  One  of  the  beams  was  shifted  about  100  Khz 
from  the  other  in  order  to  measure  the  velocity  in  the  separated  flow  region.  The 
amount  of  frequency  shift  was  decided  to  cover  the  negative  velocity  range  for  edl 
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the  test  conditions.  A  lens  of  focal  length  of  38.1cm  was  mounted  on  the  optical 
bench  to  obtain  an  intersection  volume  in  the  test  section  of  the  water  tvuinel.  The 
intensity  of  the  fringe  pattern  produced  by  the  interference  of  the  incident  beams 
was  checked  at  the  wall  furthest  from  the  test  section  using  a  microscope  objective 
located  at  the  intersection  voltime. 

FVinge  spacing  df  can  be  calculated  with  the  known  wavelength  (A)  of  the  laser 
beam  and  the  half  angle  (/c)  between  two  beams.  However,  a  direct  measurement 
was  made  using  a  rotating  disk  with  constant  angular  velocity  with  its  edge  at  the 
probe  volume  as  shown  in  Figure  12  (for  detailed  descriptions  see  Furey  and  Walker, 
1990).  With  the  known  disk  diameter  (Dr)  and  angular  velocity  (wr),  the  velocity  at 
the  edge  of  the  disk  was  calculated.  A  0.025mm  wire  attached  vertically  at  the  edge 
of  the  disk  disturbed  the  fringe  pattern  in  the  probe  volume  for  every  revolution  to 
yield  the  Doppler  frequency.  With  the  known  edge  velocity  and  Doppler  frequency, 
the  fringe  spacing,  d/,  was  calculated  as  follows, 

TT 

Uw  —  cj  ^ 

i,  =  ^  (4.3) 

Vd  ' 

where  Vx^  =  velocity  at  the  edge  of  the  rotaing  disk 
and  u}r  =  angular  velocity  of  the  disk 

The  fringe  spacing  was  measured  across  and  outside  the  probe  volume  to  check 
a  fringe  divergence,,  the  effect  of  which  was  not  significant  inside  the  probe  volume 
as  shown  in  Figure  13.  In  this  figure,  fringe  spacing  measurements  outside  the  probe 
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volume  were  included  and  marked  by  a  dotted  line.  FVinge  spacing  measured  inside 
the  probe  volume  was  approximately  5.38  /xm  with  uncertainty  of  less  than  0.8  %. 
The  uncertainty  of  each  measurement  of  the  fringe  spacing  was  accurate  up  to  0.1  % 
according  to  the  uncertainty  analysis  method  described  in  Appendix  A.  The  size  of 
the  effective  probe  volume  was  calculated  to  be  2P0/im  x  4350/xm  using  the  method 
proposed  in  Durst,  Melling,  and  Whitelaw  (1981).  Approximately  37  fringes  were 
estimated  to  exist  inside  the  probe  volume. 

Seeding  particles  were  not  necessary  because  of  many  particles  present  in  a 
water  flow. 

The  receiving  optics  of  the  same  focal  length  as  the  incident  lens  was  so  ar¬ 
ranged  that  the  scattered  beam  was  directed  to  the  photodetector.  The  focused 
beam  was  detected  by  a  UDT  PIN3D  photodiode  in  a  hybrid  cascode  circuit 
(Walker,  1984)  through  a  300  fim  pinhole. 

Details  of  the  current  LDV  setup  are  described  in  Hsu  (1989). 

4.2.S  LDV  Data  Acquisition 

The  signals  detected  by  the  photo-detector  were  processed  by  TSI  coimter-type 
processor  (Model  1990C)  with  input  conditioner,  timer,  and  analog  output  module. 
Amplitude  limit  control  was  used  to  detect  particles  producing  big  pedestals  and 
reject  the  measurement  made  on  these  particles.  The  low  limit  filter  was  used  to 
remove  the  low  frequency  component  of  the  LDV  signal  caused  by  the  Gaussian 
envelope  of  the  LDV  burst  signal  (the  pedestal  part  of  the  LDV  signed).  The  high 
limit  filter  was  used  to  eliminate  high  frequency  noise  from  the  signal  to  aid  in 
processing.  A  timer  was  used  to  measure  the  width  of  the  gate  signal  (the  period 
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corresponding  to  a  fixed  number  of  cycles  in  a  burst)  by  counting  the  number 
of  cycles  of  a  4-phase,  250  Mhz  clock  and  obtain  ±lns  accuracy  through  phase 
information.  Analog  frequency  output  was  generated  by  the  D/A  conversion  in  the 
cotmter  processor  and  sent  to  the  IBM  PC/AT  compatible  through  an  RTI-850  16 
bit  A/D  board.  By  the  above  procedure  time-averaged  signals  were  obtained  at 
each  measuring  location. 
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5.  IMAGE  DATA  PROCESSING 


5.1  Image  Capturing  System 

5.1.1  High  Speed  Video  System 

Since  the  physical  phenomenon  visualized  by  the  sequenced  multiple-wire  hy¬ 
drogen  bubble  wire  technique  was  a  complex  three-dimensional  tiurbulent  flow,  a 
high  speed  video  system  was  required  to  freeze  the  instantaneous  motion  of  the 
hydrogen  bubbles  in  a  sequence  of  instantaneous  images,  A  Kodak  Ektapro  1000 
Motion  Analyzer  was  used  to  record  the  images  on  digital  cassettes.  The  system 
consisted  of  a  main  video  processor,  two  high  intensity  cameras,  various  lenses  and 
accessories.  It  had  a  full  screen  framing  rate  of  up  to  1000  fps  (frame  per  second). 
The  frame  rate  used  in  this  research  to  capture  both  of  the  views  was  250  fps.  The 
resolution  of  image  in  a  full  screen  format  was  240  x  192  (pixels) 

5.1.2  Cameras  and  Setup 

Three  orthogonal  views  were  captured  by  two  camereis  located  at  the  side  and 
the  bottom  of  the  test  section.  The  side  view  consists  of  the  image  containing  the 
probe,  hydrogen  bubbles  and  the  wing  with  the  flow  direction  from  left  to  right.  The 
bottom  and  upstream  views  appear  on  the  same  image  using  mirror  arrangements. 
The  bottom  view  is  produced  by  looking  at  the  flow  field  from  the  bottom,  and 
shows  up  in  the  left  side  of  the  image  with  flow  direction  from  right  to  left.  The 
upstream  view  is  generated  by  the  camera  located  at  the  bottom  of  the  test  section 
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looking  at  the  flow  field  upstream  of  the  wing  through  the  mirror  installed  inside 
the  wing. 

Side  View 

Figure  10(a)  shows  two  flood  lights  located  at  the  opposite  side  of  the  test 
section  to  the  camera.  Masking  with  flat  black  plastic  electric  tape  was  made  on 
the  wall  of  the  test  section  to  allow  the  light  to  be  directed  toward  the  flow  region 
of  interest.  On  the  camera  side  of  the  wall  of  the  test  section,  masking  other  than 
the  rectangular  region  of  8.9cm  x  10.16cm  was  used  to  prevent  the  scattered  or 
reflected  light  from  various  sources  from  entering  the  camera  viewing  angle. 

The  side  view  camera  was  focused  on  the  plane  of  symmetry  and  located  far 
apart  from  the  test  section  to  avoid  the  effect  of  image  distortion.  The  field  of  view 
was  of  the  size  of  7.1cm  x  8.13cm. 

Upstream  and  Bottom  View 

In  order  to  obtain  three  orthogoniil  views  with  two  cameras,  the  upstream  and 
bottom  views  were  arranged  to  appear  on  the  same  image  plane  while  one  camera 
captured  the  side  view  in  one  image  plane.  Figure  15  shows  the  mirror  arrangement 
producing  double  images.  Three  plane  mirrors  were  installed  inside  the  mirror  box 
to  capture  the  image  from  the  flow  field.  The  approximate  size  of  7.62cm  x  12.7cm 
was  used  for  mirror  1  and  2.  These  were  inclined  at  45°  to  reflect  the  bottom  view 
image  to  the  mirror  3  of  the  size  of  7.62cm  x  25.4cm.  The  upstream  view  was 
received  by  the  same  camera  through  two  consecutive  reflections  made  by  the  45° 
slanted  mirror  inside  the  wing  and  the  mirror  3  in  the  mirror  box.  The  mirror  box 
was  located  just  below  the  test  section  to  capture  two  images  simultaneously.  The 
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average  optical  distances  for  both  of  the  views  were  designed  to  be  approximately 
equal.  The  field  of  view  of  the  upstream  view  was  approximately  7.62cm  x  7.62, 
and  that  of  the  bottom  view  was  7.62cm  x  8.9cm. 

Camera  Setup 

Two  cameras  were  employed  for  a  stereo  vision  analysis.  Ektapro  1000  system 
was  operated  in  a  full  screen  format  where  a  full  screen  was  devoted  to  the  digital 
image  from  each  imager.  The  camera  for  the  bottom  view  was  set  to  the  high 
intensity  mode,  in  which  the  sensitivity  of  the  camera  was  increased  though  the 
resolution  of  the  image  was  reduced.  Each  view  was  recorded  on  cassettes  in  an 
alternate  manner. 

After  recording  the  flow  visualization  events  on  cassettes,  the  digitized  images 
for  each  view  were  transferred  to  IBM  PC/AT  through  an  IEEEl-488  interface  and, 
finally,  stored  on  optical  disks.  A  GPIB-AT  IEEE-48S  interface  made  by  the  Na¬ 
tional  Instruments  was  employed.  The  data  transfer  program  was  developed  using 
the  C  language  for  fast  transfer.  Figure  16  shows  the  image  capturing  system. 

5.1.S  Image  Data  Formats 

The  captured  images  were  recorded  on  the  video  cassettes  in  a  digital  format. 
Each  pixel  in  an  image  has  8  bits  corresponding  to  256  gray  levels. 

Two  formats  were  used  for  the  data  transfer.  The  compressed  data  format  is 
the  one  where  one  bit  expresses  one  pixel  of  information  (2  gray  levels).  The  file 
size  for  an  image  was  5,736  bytes.  The  other  one  is  using  one  byte  per  pixel  (full 
gray  scale  image:  2®  gray  levels)  and  the  file  size  was  45,888  bytes  per  an  image. 
The  compressed  format  saved  storage  without  losing  much  image  information  by 
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thresholding  after  preprocessing  the  original  image.  This  will  be  discussed  in  later 
part  of  this  chapter.  But  the  data  transfer  time  for  a  compressed  format  image  took 
longer  than  that  for  full  gray  scale  image  since  the  preprocessing  was  performed  on 
the  digital  image  during  the  data  transfer. 

5.2  Mathematical  Background 

5.2.1  Coordinate  System 

As  shown  in  Figure  14,  the  right-handed  coordinate  system  {X,Y,Z)  and 
(U^V^W)  with  the  origin  at  the  leading  edge  of  the  wing  was  used  as  a  global 
coordinate  system.  When  processing  the  hydrogen  bubble  images,  the  image  plane 
coordinate  system  {yi,Zt)k  was  adopted  to  describe  the  individual  pixels  in  the  ib-th 
image  plane.  The  origin  of  the  image  plane  coordinate  was  the  lower  left  corner  for 
each  image  plane. 

Another  coordinate  system  employed  in  processing  the  images  was  a  homo¬ 
geneous  coordinate  system.  The  idea  of  homogeneous  coordinates  is  a  technique 
commonly  used  in  the  development  of  computer  graphics  (Riesenfeld,  1981).  The 
homogeneous  coordinates  for  the  object  points  {xo,yo,Zo)^  iii  a  globed  coordinate 
system  are  represented  by  a  4  x  1  vector  {kxoykyo^kzg^k),  where  k  is  a  nonzero 
and  arbitrary  constant.  Therefore,  there  is  no  unique  representation  of  a  physical 
point  in  this  coordinate  system.  An  image  point  is  represented  in  homogeneous 
coordinates  by  augmenting  the  2x1  physical  vector  dimension  by  one  (scale  factor) 
and  multiplying  the  physical  image  point  coordinates  by  this  nonzero  factor. 
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S.i.i  Maihematical  Model 


The  object  of  this  section  is  to  develop  a  mathematical  model  that  relates  the 
3-D  object  point  locations  measured  with  respect  to  a  global  coordinate  system 
to  the  image  point  locations  measured  with  respect  to  the  image  plane  coordinate 
system.  This  is  described  briefly  in  Figure  17(a). 

p-p  (perspective-projective)  transform 

Consider  the  pinhole  camera  model  (see  Figure  17(b)).  The  relation  between 
object  points  and  image  points  can  be  derived  using  the  constraint  that  the  image 
point,  the  center  of  projection  and  the  object  point  should  lie  on  the  same  3-D  line 
segment. 


A?(iE,  .Tc)  (^c  ^o) 


(5.1) 


Expanding  the  above, 


which  yields. 


_yo  _  fVo 
k  if  -  Xo) 
Zo  fzo 
k  if-xo) 


(5.2) 


(5.3) 


The  above  relations  are  nonlinear  in  both  of  the  object  points  and  camera 
parameter  (/  :  focal  length). 
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Using  homogeneous  coordinate  system,  the  nonlinear  relation  can  be  written 


as 


(5,4) 


where  i,- 


and 


(0  f 
P=  0  0 

V-1  0 


0  0\ 
/  0 
0  // 


Note  that  the  representation  of  P  is  not  imique. 

If  the  global  coordinate  system  is  offset  and  rotated  from  the  image  plane 
coordinate,  the  relation  becomes 


iit  =  =  Ttl,  (5.6) 

where  R*~^  is  the  transformation  matrix  which  is  related  with  rotations  and 
translations,  and  the  index  k  represents  the  expression  in  the  k-th  image  plane. 
Here,  the  final  form  of  the  transformation  matrix  T*  can  be  expressed  as  follows, 


/an 

012 

Oi3 

Ol4 

Tfc  =  1  021 

022 

023 

024 

\a3i 

032 

033 

O34 

(5.6) 


5,2.3  Camera  Calibration 

The  treinsformation  matrix  Tjt  in  Equation  (5.5)  is  not  unique  and  consists  of  12 
entities.  Each  entity  is  a  function  of  7  parameters  (3  rotational  angles,  3  translations, 
and  /)  that  determine  the  sensor  system  geometry.  With  careful  measurements  of 
the  sensor  system  geometry,  each  entity  of  the  matrix  Tjt  can  be  determined  using 
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the  expressions  derived  in  Hall  (1972).  More  practically,  this  transformation  matrix 
can  be  obtained  using  the  camera  calibration  procedure. 

A  rectangular  grid  box  with  grid  points  known  in  the  physical  coordinate  sys¬ 
tem  was  located  in  the  flow  region  of  interest  and  used  to  calibrate  the  two  cameras. 
To  find  all  of  the  eleven  parameters  (  with  setting  ch  =  1.0  lor  uniqueness),  at  least 
six  points,  which  were  non-coplanar,  were  needed  to  close  the  equations. 

Least  squares  estimation  (Q-R  Algorithm)  was  adopted  to  calculate  matrix  el¬ 
ements  out  of  more  than  twelve  equations  for  each  view.  Q-R  algorithm  is  described 
in  Appendix  B. 

Imager  calibration  coefficients  thus  obtained  are  provided  in  Table  2  for  each 
imager.  Each  cell  of  the  table  represents  a^j  in  Equation  (5.6). 

5,2.4  Stereo  Vision 

With  the  mathematical  model  described  herein,  the  transformation  from  a 
global  location  to  image  point  locations  is  always  possible  if  the  transformation 
matrix  Tk  is  known.  However,  the  inverse  transformation  is  not  unique  rmd  only 
tells  a  line  along  which  the  object  points  corresponding  to  the  image  point  should 
lie.  This  implies  that  more  than  one  calibrated  camera  are  needed  to  determine 
the  object  point  location  in  3-D  space.  Using  two  ceimeras,  the  intersection  of 
two  lines  resulting  from  each  inverse  transformation  tells  the  global  coordinate  of 
the  object  point.  Here,  the  correspondence  between  image  points  in  each  image 
plane  is  assumed.  The  correspondence  problem  is  usually  the  most  difficult  task 
in  handling  with  a  stereo  vision  application  and  will  be  discussed  separately  in  the 
later  section.  Figure  18  shows  an  example  of  a  stereo  vision  application  and  the 
constraint  between  image  points. 
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Since  Equation  (5.5)  represents  two  equations  in  a  physical  coordinate  for  each 
imager,  there  are  four  equations  and  3  unknowns  (xo,j/o,^o)  using  a  stereo  vision. 
Regrouping  Equation  (5.5)  gives  the  following  form. 


AXo  =  F 


(5.7) 


where 


(flin  -  ai3ia:ii 
ai2i  -  duiVii 
0211  -  ^2312:12 
0221  -  <l23iyi2 


0112  —  ai32®il  0113  —  01332^11 
0122  —  Oi32j/il  Oi23  —  Oi33yji 
0212  -  CL232Xi2  0-213  —  Oz^s^i2 
0222  ~  0232 yi2  0223  ~  OzzilfiZ 


and 


(Ol34a^il  ~  0114  \ 
Ol34y»l  -  0124  I 

02342^12  -  0214  I 

0234yi2  -  0224  ' 


Least  squares  estimation  needed  to  be  employed  to  obtain  the  object  pointp 


location  (a:o,yo,^o)  in  Equation  (5.7). 


5.3  Preprocessing 

In  order  to  solve  for  Equation  (5.7),  the  image  point  locations,  (a;,-,yi)jfc,  which 
correspond  to  the  object  points  in  the  global  coordinate,  should  be  determined 
first  in  each  image  plane.  In  the  present  study,  the  image  points  to  be  traced 
were  the  centroid  locations  of  the  square  shaped  area  consisting  of  small  bubbles 
generated  from  the  hydrogen  bubble  wires.  These  centroid  locations  cannot  be 
obtained  directly  from  the  image  data  recorded  in  the  cassettes  since  the  image  data 
consist  of  just  intensity  levels  (values  of  gray  scale),  f{xi,yi),  at  every  pixel  in  the 
image  plane.  Here,  /(xj,  y,)  is  the  image  function  which  represents  a  specific  image. 
Thus,  what  should  be  done  next  is  to  transform  the  intensity  level  information  into 
the  specific  locations  in  the  image  plane.  This  can  be  achieved  by  first  changing 
256  gray  levels  into  bilevel  image  which  is  black  and  white. 
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In  the  present  study,  the  preprocessing  consists  of  three  parts:  i)  the  back¬ 
ground  subtraction,  ii)  the  local  averaging  technique,  ut)  the  histogram  equaliza¬ 
tion  technique,  and  tv)  the  thresholding  technique 

'5iS.l  Background  Subtraction 

The  typical  image  function,  /o(a:i,yi)>  can  be  seen  in  Figure  19(a).  It  shows  the 
intensity  distribution  of  the  scene  with  a  wing,  a  false  wall,  hydrogen  bubbles,  and 
light  reflection  from  various  sources,  etc.  This  figure  illustrates  that  it  is  difficult  to 
find  the  reasonable  threshold  value  that  enables  the  original  image  to  be  transformed 
into  the  bilevel  image.  A  dynamic  thresholding  technique  may  be  applied  here 
to  determine  the  locally  varying  threshold  values  according  to  the  local  intensity 
values.  However,  since  the  distribution  of  the  hydrogen  bubbles  is  not  uniform  for 
the  successive  images  the  dynamic  thresholding  needs  to  be  applied  to  every  frame. 
In  practice,  it  takes  extra  amount  of  time  in  processing  the  image  with  a  dynamic 
thresholding  technique. 

Instead,  the  subtraction  of  the  nonuniform  background  from  the  original  image 
leaves  only  the  hydrogen  bubbles  in  the  field  of  view.  Ideally,  the  intensity  levels  in 
the  image  will  vary  only  due  to  the  presence  of  the  hydrogen  bubbles.  In  practice, 
the  intensity  levels  also  included  other  noise  and  the  increased  intensity  at  the  wire 
due  to  the  generation  of  the  hydrogen  bubbles. 

The  background  subtraction  can  be  written  as  follows: 


(5.8) 
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where  /*(xi,y,)  is  the  image  function  consisting  of  the  hydrogen  bubbles  only, 
foixi,yi)  is  the  original  image  function  taken  from  the  video  system,  and  /6(x,-,y,) 
is  the  image  function  for  the  background  scene. 

In  this  study,  background  images  were  captured  before  and  after  the  run  for 
every  test  with  water  flow,  and  used  to  help  process  the  data.  B’igure  19  shows  the 
effect  of  the  background  subtraction.  In  the  original  image  as  shown  in  Figure  19(a), 
the  intensity  representing  the  hydrogen  bubbles  seemed  buried  in  other  higher  level 
of  intensity  mainly  due  to  the  background.  However,  the  local  peaks  of  intensity 
showing  the  existence  of  the  hydrogen  bubbles  appear  after  background  subtraction 
as  shown  in  Figure  19(c). 

5.S,£  LocaLAveraging  Technique 

The  image  function,  A(a:i,yi),  obtained  after  subtracting  the  background  sub¬ 
traction,  still  contained  noise  imposed  on  the  intensity  levels,  which  needed  to  be 
eliminated.  One  source  that  arises  from  an  electronic  sensor  appears  as  “snow”, 
which  is  random  and  uncorrelated.  Since  this  type  of  noise  shows  extreme  change 
in  intensity  between  the  adjacent  pixels,  a  local  averaging  technique  was  used  to 
remove  it. 

In  this  technique,  the  intensity  value  that  shows  the  extreme  difference  is  re¬ 
placed  by  a  local  average  calculated  using  neighboring  pixels.  This  nonlinear  oper¬ 
ation  can  be  described  by 


y(x,-,yi)  =  I  fhi^i^yi)\>'^  (5,9) 

Xhi^i^yi)  otherwise, 


where  T  is  a  threshold  value  dependent  on  the  noise  variance. 
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The  present  study  used  8  neighboring  pixels  for  the  local  averaging.  Figure  20 
describes  the  windows  used  for  calculating  the  local  average  in  the  vicinity  of  a 
pixel.  The  center  point  in  the  3x3  window  was  excluded  to  eliminate  extreme 
values  in  the  averaging.  If  too  many  pixels  are  subject  to  local  averaging,  the  image 
is  further  smoothed  out  with  the  adverse  effect  of  smearing  the  image. 

S.S.S  Histogram  Equalization  Technique 

For  the  image  function,  g(xi,yi)^  the  intensity  levels  of  the  image  very  near  the 
leading  edge  of  the  wing  are  low  when  compared  to  those  at  other  pixel  locations. 
The  intensity  of  the  reflected  light  from  the  wing  has  almost  the  same  level  as 
those  of  the  hydrogen  bubbles.  We  lose  contrast  in  this  region,  and  therefore,  need 
enhance  the  image  by  further  processing. 

As  an  contrast  enhancing  technique,  a  histogram  equalization  technique  was 
adopted  since  a  uniform  distribution  of  gray  levels  makes  equal  use  of  all  the  quan¬ 
tization  levels  and  enhances  low-contrast  information.  The  algorithm  for  this  tech¬ 
nique  is  presented  in  Figure  21,  and  will  be  described  as  follows. 

If  P  is  the  total  number  of  pixels  in  an  image,  and  N  is  the  number  of  avail¬ 
able  gray  levels,  one  must  have  P/N  pixels  at  each  gray  level  for  a  perfectly  flat 
histogram.  Therefore,  for  the  intensity  level  (G)  with  the  number  of  occurrences  k 
times  greater  than  the  average  level,  that  particular  level  should  be  redistributed 
to  k  different  levels,  i.e.  (Gi,Gjt). 

There  are  several  possible  ways  to  redistribute  the  pixel  intensity  (Pavlidis, 
1982).  But,  in  the  present  study,  an  empirical  distribution  function  was  computed 
by  adding  the  relative  frequency  values  and  used  for  reassigning  a  new  intensity 
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level  at  low  contrast  re^on.  The  distribution  curve  illustrates  that  the  low  contrast 
region  may  be  enhanced  by  stretching  the  low  intensity  range  so  that  these  values 
can  be  reassigned  to  higher  levels  of  intensity  according  to  the  distribution  curve. 

The  steps  presented  below  are  the  procedure  used  for  the  histogram  equaliza¬ 
tion: 

(a)  Compute  the  histogram  iH(G))  of  the  image  gray  level  valves, 

(b)  Obtain  a  distribution  curve  by  accumulating  the  histogram  values, 

(c)  Use  this  distribution  curve  for  the  gray  level  transformation 

h{xi,yi)  =  D{fhixi,yi))  (5.10) 

where  D  is  the  operator  performing  the  redistribution  of  the  intensity  level. 
5.3.4  Thresholding  Technique 

The  final  step  to  obtain  a  bilevel  image  employs  a  thresholding  technique.  In 
this  technique,  the  intensity  value  at  each  pixel  is  compared  to  a  given  threshold 
value,  and  the  pixel  is  assigned  to  one  of  two  categories  depending  on  whether  the 
threshold  value  is  exceeded  or  not. 

The  threshold  value,  in  general,  is  selected  from  the  intensity  histogram.  For 
example,  if  the  histogram  has  two  peaks,  each  representing  the  bright  area  and  the 
dark  area  in  the  image,  a  reasonable  choice  will  be  a  point  between  the  two  peaks. 
However,  the  histogram  of  the  intensity  does  not  always  show  two  peaks,  and  the 
intensity  level  for  each  image  can  vary  due  to  the  limitations  of  the  Ektapro  1000 
System.  Therefore,  we  could  not  obtain  a  uniform  value  for  the  threshold  for  the 
image  used  in  data  reduction. 
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Instead,  the  threshold  level  was  established  automatically  in  the  data  reduction 
program  for  the  image,  based  on  a  table  for  the  percentage  of  “ON”  pixels  for  a 
specific  scene.  In  order  to  do  this,  the  image  sequences  were  recognized  using  the 
location  of  the  hydrogen  bubbles  in  the  flow  field.  This  was  possible  because  the 
five  wires  generate  hydrogen  bubbles  at  different  times  and  locations. 

This  technique  can  be  briefly  described  by 


if  /»(x,-,  yi)  >  Th 
otherwise. 


(5.11) 


where  y,-)  is  the  image  function  of  bilevel  intensity,  and  Th  is  a  selected  thresh¬ 
old  value.  Figure  22  shows  examples  of  bilevel  images  after  thresholding  analysis. 


IMAGE  DATA  PROCESSING 


42 


5.4  Motion  Analysis 


5.4>1  Introduction 

The  bilevel  images  constructed  after  preprocessing  are  made  up  of  the  pixels 
with  two  gray  levels  (“ON”  and  “OFF”  pixels).  Hydrogen  bubbles  Me  represented 
by  the  “ON”  pixels.  Figure  22  illustrates  the  bilevel  image. 

The  three-dimensional  analysis  of  the  motion  of  the  flow  used  in  this  study 
determines  the  three-dimensional  instantaneous  velocity  field.  To  do  this,  the  gray 
level  information  in  the  images  should  be  transformed  into  the  information  of  bubble 
location.  The  displacement  of  a  particular  bubble  region  for  successive  frames  is 
used  to  calculate  the  corresponding  instantaneous  velocity. 

This  can  be  obtained  by  applying  the  following  algorithms  successively. 

1)  Centroid  search  routine 

In  this  routine,  a  bubble  element  numbering  is  done  by  tracing  “ON”  pixels 
in  the  j-th  frame.  For  the  i-th  bubble  element,  the  location  of  the  centroid 
is  determined. 

2)  Correspondence  Algorithm 

This  algorithm  determines  the  correspondence  of  the  image  points  between 
the  image  planes.  The  same  bubble  element  number  is  assigned  to  the 
corresponding  bubble  elements. 

3)  Matching  algorithm 

In  this  routine,  matched  pairs  of  the  same  bubble  elements  between  two 
successive  frames  (the  ;-th  and  (j+l)-th.  frames)  are  determined.  Two 
criteria  are  applied  to  determine  the  matched  pairs. 
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4)  Split-merging  algorithm 

In  this  algorithm,  an  attempt  is  made  to  retrieve  the  bubble  elements  lost 
due  to  the  merging  of  them  in  the  image  plane. 

Using  the  above  algorithms,  the  corresponding  »-th  element  in  both  image 
planes  is  found  and  the  element  in  the  j-th  frame  which  is  matched  to  the  i-th 
bubble  element  in  the  (j+1  j-th  frame  is  determined.  Thus,  the  physical  location  of 
the  t-th  bubble  element  can  be  calculated  using  Equation  (5.7)  with  the  constraints 
on  the  correspondence  between  the  bubbles.  The  difference  in  the  centroid  locations 
of  the  i-th  element  between  the  ;-th  and  (ji-l)-ih  frames  can  be  used  to  calculate 
the  displacement  vectors  of  the  individual  bubble  elements  between  frames. 

5,4-2  Ctntroid  Search  Routine 

As  shown  in  Figure  22,  we  needed  some  criteria  to  determine  the  bubble  el¬ 
ements  and  numbering.  The  procedure  to  search  the  centroid  location  of  edl  the 
bubble  elements  is  as  follows. 

The  groups  of  the  “ON”  pixels  were  searched,  row  by  row,  starting  from  the 
origin  in  the  image  plane  to  the  last  pixel.  Once  the  first  “ON”  pixel  was  found, 
the  neighboring  pixels  were  checked  if  there  were  other  “ON”  pixels,  and  this  was 
applied  to  all  the  neighboring  pixels  until  there  were  no  more  “ON”  pixels.  The 
noise  was  still  imposed  on  the  bilevel  image  with  all  the  preprocessing  techniques 
described  in  the  previous  section  and  the  bubble  elements  with  smaller  numbers  of 
pixels  appeared  irregularly  from  frame  to  frame.  Therefore,  the  groups  with  very 
small  numbers  of  the  “ON”  pixels  were  not  considered  to  be  bubble  elements.  This 
caused  some  bubble  elements  to  be  lost  during  the  centroid  search,  but  resulted 
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in  less  error  in  the  computer  processing  than  the  inconsistent  appearance  of  some 
images  from  frame  to  frame.  After  the  first  bubble  element  was  found,  the  routine 
went  to  the  previous  row  where  the  first  “ON”  pixel  was  found,  and  restarted  to 
search  the  next  “ON”  pixel  in  the  same  row.  This  routine  was  applied  to  the  end  of 
the  row  to  seairch  all  the  bubble  elements  that  had  a  pixel  or  pixels  belonging  to  the 
same  row.  The  bubble  elements  found  in  sequence  were  assigned  unique  numbers 
in  an  increasing  order. 

In  searching  the  bubble  elements,  the  “ON”  pixels  belonging  to  a  bubble  el¬ 
ement  may  be  well  closed  and  do  not  show  any  confusion  with  any  other  bubble 
elements.  However,  this  did  not  happen  always.  Sometunes,  two  bubble  elements 
merged  to  one,  while  one  element  split  into  two  or  more  bubble  elements  due  to 
strong  vortex  motion  or  the  diffusion  of  the  hydrogen  bubbles.  In  addition,  the 
bubble  elements  at  different  locations  in  the  physical  coordinate  could  appear  in 
the  same  location  in  the  image  plane.  Therefore,  the  criteria  on  the  connectivity  of 
the  pixels  to  form  a  bubble  element  should  not  be  rigid  to  take  the  above  cases  into 
account.  Figure  23  shows  the  example  of  the  bubble  element  identification.  Here, 
the  pixels  connected  diagonally  to  another  group  of  pixels  were  treated  eis  separate 
bubble  elements  as  shown  in  Figure  23(b).  In  this  figure,  it  is  illustrated  that  some 
isolated  “ON”  pixels  connected  diagonally  to  a  group  of  “ON”  pixels  axe  treated  as 
not  belonging  to  that  group. 

From  the  next  row  to  the  last,  the  algorithm  checked  if  the  “ON”  pixels  be¬ 
longed  to  the  previous  bubble  elements  and,  if  so,  were  excluded.  The  same  routine 
was  applied  to  all.  the  rows  in  each  image  plane  until  all  the  bubble  elements  could 
be  identified. 
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For  eadi  bubble  element  found  in  this  routine,  the  centroid  location  could  be 


calculated  using  the  pixel  locations  in  the  image  plane  as  follows, 

Xci  = 


Vci  = 


Ni 

ENi 

kUVk 

Ni 


(5.12) 


where  Ni  is  the  number  of  pixels  of  the  i-ih  bubble  element 


5.4>2  Correspondence  Algorithm 

After  the  centroid  locations  of  the  bubble  elements  existing  in  the  image  plane 
were  found,  it  remained  to  determine  the  coordinates  of  the  centroid  locations. 
To  obtain  the  locations  of  the  bubble  elements  in  the  three-dimensional  space, 
the  correspondence  of  the  bubble  elements  between  the  image  planes  should  be 
recognized.  This  is  the  most  difficult  problem  in  processing  the  stereo  vision  images. 

As  expressed  in  Equation  (5.7),  knowing  the  corresponding  points  in  each  image 
plane,  the  object  point  in  the  global  coordinate  system  could  be  determined  by 
applying  least  squares  analysis  to  these  equations.  However,  not  all  the  locations 
are  required  to  estimate  Xo^  The  object  point  Xo  is  determined  using  three  out  of 
four  equations,  assuming  that  the  three  equations  were  linearly  independent. 

In  a  homogeneous  coordinate  system.  Equation  (5.7)  becomes 


CXo  =  0 


(5.13) 


where  C  = 


/  —  Oni 

o-mVii  -  0121 
02313;, 2  -  0211 
\a231*/i2  -  0221 


01323;,1  —  Oii2 
Ol32yil  -  0122 
02323:,2  —  0212 
0232yi2  -  0222 


Ol332;il  —  Oii3 
Ol33yil  -  0123 
02333;,2  -  0213 
0233yi2  ~  0223 


01343;,1  —  Oii4  \ 
Ol34y»l  -  0124 
023437, ‘2  -  0214 
0234  y«2  ~  0224  / 
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Equation  (5.13)  requires  that  every  row  (or  column)  of  is  orthogonal  to  Xo- 
The  algebraic  significance  of  Equation  (5.13)  is  that  the  column  must  not  span 
which  indicates  that  the  determinant  of  matrix  C  should  be  zero.  Therefore, 

det  \C\  =  0  (5.14) 

which  yields  the  following  equation. 


miXiiXi2  +  m2yiiy<2  +  JT^s^.iyiz  + 

(5.15) 

+m5a;ii  +  me®, -2  +  m7y,i  +  m%yi2  +  mg  =  0. 

The  coefficients  mj’s  given  in  Equation  (5.15)  are  functions  of  the  sensor 
calibration  parameters  and  derived  in  Appendix  C.  The  geometric  significance  of 
Equation  (5.15)  is  that,  given  a  point  in  one  image  plane  e.g.,  (a:i,,y«i)>  the  set  of 
corresponding  image  points,  {xi,yyij),  in  the  second  image  plane  is  constrained  to 
lie  along  a  line  formed  by  Equation  (5.15).  The  constraint  on  the  correspondence 
problem  is  illustrated  in  Figure  18  and  Figure  24. 

The  constraints  given  in  Equation  (5.15)  cannot  solve  the  correspondence  prob¬ 
lem  uniquely  because  any  points  near  that  line  are  the  possible  candidates  for  the 
corresponding  point.  In  our  study,  the  sequence  of  the  generation  of  bubble  from 
the  wires  was  used  for  another  constraint  to  determine  the  image  point  in  one  image 
plane  corresponding  to  the  point  in  other  image  plane. 

At  the  instant  of  pulsing  to  the  wire,  the  bubble  elements  were  generated 
simultaneously  at  some  band  of  Z  coordinate  in  the  image  plane  representing  a 
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upstream  view  and  there  were  five  bands  corresponding  to  each  wire.  The  bubble 
elements  fotmd  during  the  centroid  search  routine  described  in  the  previous  section 
would  be  pven  their  unique  numbers  from  the  bottom  to  the  top  in  the  side  view 
image  plane.  The  corresponding  bubble  elements  were  easily  found  in  the  upstream 
view  image  plane  because  the  bubble  elements  should  be  located  at  the  first  Z 
band  and  the  bubble  elements  which  satisfied  Equation  (5.15)  should  lie  near  the 
line  which  is  almost  normal  to  that  band  of  constant  Z  line.  This  means  that 
the  corresponding  bubble  elements  should  be  located  near  the  intersection  of  a 
corresponding  line  and  a  band  of  constant  Z.  After  all  the  corresponding  bubble 
elements  were  found  in  the  j-th  frame,  the  corresponding  bubble  elements  were 
always  given  the  same  number.  From  (j-hl)-th.  frame  on,  this  bubble  number  was 
used  to  trace  the  corresponding  bubble  elements. 

Since  we  know  the  sequence  of  bubble  generation,  the  next  expected  Z  band 
in  the  upstream  view  image  plane  can  be  searched  for  the  correspondence  at  the 
bubble  generation  from  the  next  wire,  and  the  numbers  that  followed  the  previous 
numbering  would  be  distributed  to  the  newly  found  bubble  elements  until  the  whole 
cycle  was  completed. 

5. 4-3  Matching  Algorithm 

We  applied  the  matching  algorithm  to  the  successive  frames  (j-th  and  (j+1)- 
th  frames)  for  the  same  image  plane  to  find  the  same  bubble  elements  .as  time 
proceeded.  The  criteria  on  the  matched  bubble  elements  were  using  the  area  of  the 
individual  bubble  elements  in  the  (j+l)-ih  frame,  which  were  expected  to  be  at 
some  pixel  distance  from  the  centroid  location  in  the  ;-th  frame. 


IMAGE  DATA  PROCESSING 


48 


The  distance  between  the  matched  bubble  elements  in  the  j-th  and  (j+1  )-th 
was  restricted  by  the  estimated  maximum  flow  velocity.  Thus,  only  the  bubble 
elements^  the  centroid  locations  of  which  were  within  certain  range  of  the  distance 
from  the  previous  centroid  location,  were  investigated  to  search  the  matched  bubble 
element. 

The  maximum  expected  pixel  distance  was  approximately  3  pixels,  which  was 
obtained  by  considering  the  screen  resolution,  the  field  of  view  in  the  absolute 
scale,  and  the  maximum  expected  flow  speed.  The  calculation  of  the  pixel  distance 
is  presented  in  Table  4. 

5.4-4’  Split-Merging  Detection  Algorithm 

All  the  bubble  elements,  which  at  first  were  independent  of  each  other,  could 
merge  together  as  the  bubbles  approached  the  wing.  Especially,  the  bubble  elements 
in  the  strong  vortex  region  were  easily  stretched  and  mixed  with  th?  others,  and, 
in  the  long  run,  became  not  discernible  with  the  others.  Once  the  bubble  elements 
entered  the  strong  vortex  region,  it  was  difficult  to  retrieve  those  bubble  elements  as 
separate  ones.  In  addition,  though  the  bubble  elements  were  not  actually  merged, 
it  could  happen  that  several  bubble  elements  in  3-D  space  appeared  to  be  just  one 
bubble  element  in  2-D  image  plane. 

Therefore,  efforts  were  made  to  extract  more  information  in  the  region  very 
near  the  wing  by  retrieving  the  bubbles  lost  due  to  merging  and  other  reasons.  In 
particular,  as  illustrated  in  Figure  24,  for  the  bubble  elements  that  took  apart  in 
one  image  plane,  but  appeared  merged  in  other  image  plane,  we  could  retrieve  the 
lost  information  of  the  centroid  though  it  was  subject  to  have  some  uncertainties. 
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The  algorithm  for  the  detection  of  the  split-merging  of  the  bubble  elements  is 
provided  as  a  flowchart  in  Figure  25. 

5.4-5  Instantaneous  Velocity  Field 

After  a  cycle  of  the  motion  analysis  algorithms  described  above  is  completed,  a 
set  of  bubble  element  locations  in  the  global  coordinate  is  acquired.  As  the  process 
goes  on,  the  trajectories  of  each  bubble  element  in  the  global  coordinate  are  obtained 
as  a  solution  set. 

The  velocity  vector  in  the  (j-fl )-th.  frame,  Uj+i,  can  be  written  as  follows, 

(££j+i  -  £»;•) 

-  ~~Tf 

where  t/  is  the  sampling  interval  between  j-th  and  (j+l)-tlci  frames. 

Uncertainty  Consideration 

As  provided  in  Table  4,  the  maximum  expected  pixel  distance  in  the  axial 
direction  for  the  given  test  condition  was  3  pixels  between  successive  frames.  This 
indicates  that  the  uncertainty  in  the  axial  velocity  is  at  least  16.7%  since  the  centroid 
locations  of  each  bubble  element  are  always  subject  to  the  error  of  half  pixel  distance. 

One  way  to  reduce  such  a  large  uncertainty  in  velocity  calculation  is  to  use  sev¬ 
eral  successive  frames  and  compute  the  averaged  velocity  of  certain  bubble  element 
in  these  frames,  and  then  use  this  average  value  for  the  velocity  in  certain  frame. 
Using  five  frames  to  obtain  the  averaged  velocity,  the  uncertainty  was  reduced  to 
7.5%. 
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6.  EXPERIMENTAL  RESULTS 


6.1  Introduction 

This  chapter  presents  the  experimental  results  of  the  hydrogen  bubble  visu¬ 
alization,  LDV  measurements,  and  three-dimensional  PIDV  (Particle  Image  Dis¬ 
placement  Velocimetry)  using  a  high  speed  video  system.  Experimental  results  are 
provided  using  a  right-handed  coordinate  system  for  physical  plane  data,  and  an 
image  coordinate  system  for  image  plane  data  as  described  in  Figure  14.  The  un¬ 
certainties  of  the  experimental  results  were  estimated  by  the  method  discussed  in 
Appendix  A,  and  are  tabulated  in  Table  2. 

The  first  part  of  this  chapter  consists  of  a  general  description  of  the  junction 
flow  using  hydrogen  bubble  flow  visudization.  Several  sequences  of  hydrogen  bubble 
pictures  are  presented  for  illustration. 

LDV  measurements  are  provided  to  show  the  general  mean  flow  characteristics 
in  the  plane  of  symmetry  upstream  of  the  wing.  These  are  compared  with  the  wind 
tunnel  data  of  Devenport  and  Simpson  (1987). 

Next,  the  three-dimensional  instantaneous  velocity  vector  field  in  the  junction 
region  is  presented,  which  was  obtained  by  the  3-D  PIDV  method. 


EXPERIMENTAL  RESULTS 


51 


6.2  Test  Conditions 

The  requirements  for  the  test  conditions  for  the  present  study  were  that  the 
potential  flow  be  steady  and  uniform,  and  the  approach  boundary  layer  should  be 
turbulent. 

For  flow  visualization  using  the  hydrogen  bubble  technique,  the  flow  speed  was 
restricted  due  to  the  production  of  cavitation-induced  air  bubbles,  which,  at  higher 
speed,  interfered  with  the  hydrogen  bubbles.  For  best  flow  visualization  results,  the 
flow  speed  was  approximately  9.1cm/5ec  at  XjT  =  -5.33,  which  corresponded  to 
Rtt  —  330  at  XJT  —  -2.15.  Here,  T  is  the  maximum  thickness  of  the  wing  and 
Re$  is  the  Reynolds  number  based  on  the  momentum  thickness  of  the  oncoming 
boundary  layer.  However,  since  at  this  flow  speed  the  boimdary  layer  was  barely 
turbulent,  a  higher  speed  was  also  used  with  LDV  measurements  for  comparison. 
For  the  higher  speed.  Ret  =  1100  was  obtained  at  XjT  =  -2.15.  Test  conditions 
are  listed  in  Table  1. 

6.3  Flow  Visualization  Results 

Flow  visualization  using  a  hydrogen  bubble  wire  technique  was  used  in  the  nose 
region  of  the  wing-body  junction.  Tests  were  conducted  at  low  water  speed  {Re^  = 
330),  and  flow  visualization  results  were  taped  using  a  standard  VHS  camcorder. 
Photographs  presented  in  this  section  were  obtained  via  a  Nikon  2020  35mm  camera 
either  directly  from  the  flow  or  from  a  TV  screen  displaying  a  frozen  frame  of 
a  recorded  VHS  tape  on  a  Panasonic  Model  AG1960  super  VHS  VCR.  An  FDL 
(Fluorescent  to  Daylight)  filter  was  used  for  the  TV  photographs  for  improved 
color  balance. 
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Figure  26  shows  a  sequence  in  which  inverted  time-lines  generated  every  second 
were  used  as  flow  tracers.  The  location  of  the  hydrogen  bubble  wire  was  at  6.35cm 
upstream  of  the  leading  edge  of  the  wing  {XjT  =  —0.443).  It  should  be  noted  that 
the  wire  was  not  exactly  vertical  to  the  floor  (about  6“  inclined  toward  the  wing), 
but  still  in  the  plane  of  symmetry.  The  dark  area  near  the  wing  is  in  the  shadow  of 
the  wing.  The  field  of  view  (from  the  wire  to  the  leading  edge  of  the  wing)  is  about 
6.35cm  X  7.20cm.  The  disturbed  flow  near  the  top  of  the  view  was  caused  by  the 
interference  of  the  wire  junction  of  the  probe  with  the  flow. 

Figure  27  shows  a  sequence  of  pictures  which  captured  the  flow  near  the  saddle 
point  formed  between  two  vortices.  The  location  of  the  wire  was  at  5.50cm  upstream 
of  the  leading  edge  of  the  wing.  The  field  of  view  is  about  5.50cm  x  7.20cm.  Near 
the  top,  the  wire  had  a  foreign  material  such  as  dust  or  small  bubbles,  so  that  no 
bubbles  were  generated  on  that  spot. 

Figure  28  shows  the  time-line  patterns  of  the  various  flow  structures  found 
in  this  region.  The  wire  was  at  7.62cm  upstream  of  the  leading  edge  of  the  wing 
{X/T  =  —0.531).  The  field  of  view  is  7.62cm  x  7.94cm. 

6.4  Results  of  LDV  Measurements 

LDV  measurements  were  made  in  the  plane  of  symmetry  upstream  of  the  wing. 
Figure  29  shows  the  measurement  locations. 
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6.4.1  Characterisiics  of  Approaching  Boundary  Layer 

The  characteristics  of  the  approaching  boundary  layer  were  investigated  at 
X/T  =  -2.15  using  LDV  measurements  (Figure  11  shows  the  setup),  and  the  mean 
flow  data  at  this  location  were  used  as  parameters  for  comparison  with  the  wind 
timnel  data  of  Devenport  and  Simpson  (1987).  The  boundary  layer  characteristics 
at  this  location  are  tabulated  in  Table  1. 

The  momentum-thickness  Reynolds  numbers  {Reo)  measured  at  this  location 
were  330  and  1100,  respectively,  under  the  nominal  test  conditions.  These  corre¬ 
spond  to  the  boundary  layer  thicknesses  (^99.5)  of  40.6mm  =  0.283)  for  the 

Res  =  330  case,  and  27.3mm  (^99.5/T  =  0.190)  for  the  Res  =  1100  case.  The  mean 
velocities  measured  at  X/T  -  -5.312  were  used  as  reference  velocities  {Urtf)i 
Urtf  of  9.07cm/s  and  44.75cm/s  were  obtained,  respectively. 

Figu”';  30  shows  the  comparison  of  the  law  of  the  wall  plots  measured  ut  X/T  — 
—2.15  for  each  test  condition  including  Devenport’s  wind  tunnel  data.  Three  data 
sets  coalesced  fairly  well  in  the  logarithmic  region.  The  Reynolds  number  effect 
shows  up  in  the  outer  region  of  this  plot  where  the  highest  Reynolds  number  case 
deviates  up  from  the  slope  of  the  logarithmic  region.  The  present  data  do  not  have 
enough  data  points  in  this  region,  but  have  the  trend  showing  the  Reynolds  number 
effect. 

Figure  31  shows  the  distribution  of  the  axial  turbulent  stresses  across  the 
boundary  layer  for  both  of  the  test  conditions.  It  should  be  noted  that  the  tur¬ 
bulence  level  in  the  external  flow  was  high  as  shown  in  the  hot-film  measurement 
data  in  Figure  3. 


EXPERIMENTAL  RESULTS 


54 


Figures  32  -  37  show  the  procedures  of  obtaining  the  local  skin  friction  coef¬ 
ficients  (Cf)  at  X/T  =  -2,15.  For  the  high  Reynolds  number  case,  the  local  skin 
friction  coefficients  of  0.0053  and  0,0078  were  obtained  using  the  “cross-plot  of  log 
region”  method  and  the  “wall-profile”  method,  respectively.  Law  of  the  wall  plots 
were  made  for  both  of  the  test  conditions. 

6.4^B  Mean  Flow  Data 

The  mean  velocity  profiles  at  all  the  measurement  locations  are  plotted  in 
Figure  38  and  Figure  39  for  each  test  condition.  The  mean  flow  profile  at  X/T  = 
—0.86  was  still  very  similar  to  that  of  the  two-dimensional  turbulent  boundary  layer. 
As  the  flow  approached  the  wing,  the  effect  of  the  primary  junction  vortex  appeared 
in  the  outer  part  of  the  boundary  layer  for  both  of  Res  =  330  and  Reg  =  1100  cases. 
The  flow  in  the  outer  layer  was  observed  to  be  accelerated  faster  than  the  external 
flow  due  to  the  action  of  the  strong  junction  vortex.  This  happened  at  all  the 
measuring  points  from  X/T  =  —0.67  to  the  point  nearest  to  the  leading  edge  of  the 
wing  for  Reg  =  1100,  and  from  X/T  =  —0.40  on  for  Reg  =  330.  This  indicates  that 
the  influence  of  the  junction  vortex  reaches  far  upstream  at  higher  speed. 

The  backflow  reached  a  maximum  of  —OAlUref  at  X/T  =  -0.125  and  Y/T  = 
0.0078  for  the  Reg  —  330  case  and  of  -0.48C^re/  at  X/T  =  -0.075  and  Y/T  = 
0,0032  for  the  Reg  =  1100  case.  Comparing  with  Devenport’s  wind  tunnel  data  of 
—QAWrtf  at  X/T  =  —0.153  and  Y/T  =  0.0047,  the  present  data  shows  that  the 
maximum  backflow  scaled  with  Uref  <hd  not  change  much  at  higher  Reg. 

Figure  40  shows  the  boundary  layer  profiles  at  all  the  measuring  points  in  the 
plane  of  symmetry  upstream  of  the  wing.  This  plot  was  constructed  using  only  U- 
component  velocities.  The  region  surrounded  by  a  solid  line  indicates  the  bimodal 
zone  found. 
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In  Figure  41,  three  data  sets  with  different  Ree  were  plotted  together  to  teike  the 
effect  of  the  Reynolds  number  into  accoimt.  In  these  plots,  Y-location  was  scaled 
with  the  approach  boundary  layer  thickness  6,  while  X-locations  of  measurement 
were  scaled  with  T. 

6.4’S  Some  Turbulence  Data 

The  variations  across  the  boundary  layer  were  obtained  at  all  the  mea¬ 
surement  locations  and  plotted  with  Devenport’s  wind  tunnel  data  as  shown  in 
Figure  42.  These  data  show  how  the  flow  fluctuations  near  the  junction  vary  as 
the  wing  is  approached.  At  all  the  measuring  points,  the  characteristics  of  the  flow 
deviated  from  the  typical  turbulent  stress  distribution  in  the  boundary  layer  due  to 
the  presence  of  the  wing. 

Here,  the  U-component  fluctuations  (  )  increased  as  the  flow  ap¬ 

proached  the  wing.  The  magnitude  of  fluctuations  reached  a  maximum  approxi¬ 
mately  at  X/T  =  —0.125,  or  —0.150,  and  then  reduced  to  very  small  value  in  the 
vicinity  of  the  wing.  In  these  plots,  the  variations  show  similar  behavior  across 
the  boundary  layers  for  all  the  test  cases.  In  particular,  in  the  case  of  Reg  =  1100 
and  Devenport’s  data,  where  both  of  them  were  comparatively  high  Reynolds  num¬ 
ber  cases,  the  local  maximum  in  at  high  Y  location  was  observed  first, 

and  then  another  peak  in  at  low  Y  location  followed  as  the  wing  was  ap¬ 

proached.  Figure  43(a)  and  Figure  43(b)  show  the  contour  plots  of  for  Reg  =  330 
and  Reg  =  1100,  respectively. 

The  velocity  histograms  across  the  boundary  layer  at  each  measuring  station 
where  vortex  motion  was  considered  to  be  approximately  centered,  are  presented  in 


EXPERIMENTAL  RESULTS 


56 


Figure  44.  The  tic  marks  in  these  plots  show  zero  velocity  points  in  the  histogram. 
Here,  one  can  see  the  bimodal  (double-peaked)  velocity  histograms  in  some  XjT 
and  YfT  range.  Note  that  these  bimodal  velocity  histograms  appeared  in  both  of 
the  test  cases. 

Higher-order  turbulent  results 

The  skewness  factor  (Su)  and  flatness  factor  (Fu)  for  u  were  calculated  from 
the  velocity  histograms  P{U)  at  all  the  measurement  locations. 

The  skewness  factors  {Su)  and  flatness  factors  (Ft)  are  given  by 
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where  the  third  and  fourth  moments  can  be  obtained  by  setting  n  =  3  and  4  from 
the  following  expression, 
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The  results  are  presented  in  Figure  45  -  Figure  48  for  each  test  condition. 


6,^-4  Switching  Frequency  Analysis  (SFA) 

As  pointed  out  by  Devenport  and  Simpson  (1987, 1990),  the  flow  in  this  region 
appears  to  have  two  distinct  modes:  a  backflow  mode  and  a  zero  flow  mode.  In 
the  backflow  mode,  the  flow  in  the  separation  region  reaches  far  upstream  forming 
a  large  vortical  structure,  while  the  fluid  in  the  zero  flow  mode  flows  downstream 
reducing  the  extent  of  the  separated  region. 
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For  the  LDV  data  which  had  bimodal  velocity  histograms,  the  average  fre¬ 
quency  of  switching  (/,„,)  from  one  mode  of  flow  to  another  was  calculated  to  study 
the  relations  among  the  experimental  parameters,  especially  between  the  switching 
frequency  and  the  flow  structure. 

A  Single  Threshold  Level  SFA  with  6Tm;„ 

This  analysis  method  employed  a  single  threshold  between  two  peaks  in  the 
velocity  histogram  and  a  minimum  switching  time  (STmin)  as  parameters  to  deter¬ 
mine  the  switching  of  flow  mode.  In  this  technique,  flow  mode  was  considered  to  be 
switched  when  time  series  of  LDV  signals  encountered  a  threshold  level  and  lasted 
more  than  ST^in  in  one  flow  mode. 

Figure  49  illustrates  some  results  obtained  by  this  analysis  method.  Eleven 
threshold  levels  were  selected  between  two  peaks  and  STmin  also  varied  from  0.075 
sec.  to  0.20  sec.  As  seen  in  these  plots,  the  average  switching  time  (T,u))  depends 
on  the  two  selected  parameters.  However,  it  is  notable  that  Tsw  did  not  change 
much  with  a  selected  threshold  level  when  two  flow  modes  existed  2ilmost  at  the 
same  rate. 

Three  Threshold  Level  SFA  with  Digital  Filtering 

In  this  method,  the  time  series  of  velocity  data  were  filtered  digitally  to  remove 
the  high  frequency  turbulent  fluctuations  which  were  considered  not  to  be  associated 
with  the  change  of  the  flow  mode.  Four  different  kinds  of  FIR  (Finite  Impulse 
Response)  filters  were  designed  according  to  N  (number  of  filter  taps)  and  3dB 
points,  arid  compared  to  each  other.  The  characteristics  of  each  digital  filter  are 
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presented  in  Figure  50.  Figure  51  illustrates  how  each  filter  works  on  LDV  signals 
in  the  bimodal  region. 

After  passing  through  the  digital  filter,  the  velocity  histogram  was  obtained  and 
smoothed  further  to  calculate  the  extremal  locations  in  the  histogram  without  the 
influence  of  the  local  variation  of  the  shape  of  the  histogram.  This  was  confirmed  at 
every  data  point  by  the  visual  inspection  through  the  program  developed  to  produce 
an  interactive  graphics  during  computation.  Two  local  maximum  locations  (vmoxi 
and  Vfnaxi)  minimum  location  (umm)  were  obtained  in  a  velocity 

histogram  which  had  a  bimodal  shr^pe. 

Three  threshold  levels  between  minimum  and  maximum  velocity  were  selected 
to  trigger  the  switching  between  the  modes.  Since  the  shape  parameters  of  a  bi¬ 
modal  histogram,  such  as  the  heights  of  two  peaks  and  the  widths  between  the  local 
extremes,  are  dependable  factors,  the  threshold  levels  were  selected  to  include  the 
effect  of  those  parameters  as  below.  1 
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where  hi  and  hz  are  the  heights  of  each  peak  in  the  velocity  histograms.  Figure  52  is 
the  schematic  drawing  describing  the  parameters  employed  in  Equation  (6.2).  The 
constants  Ci  and  C3  in  Equation  (6.2)  were  selected  to  give  better  visual  results  on 
the  values  of  and  Vt^ ,  and  the  values  of  0.5  for  each  constant  worked  well  after  a 
series  of  numerical  experiments.  The  typical  example  of  the  triggered  LDV  signals 
using  three  threshold  values  is  given  in  Figure  53. 
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Figure  54  shows  the  variations  of  the  probability  density  function  of  the  switch¬ 
ing  frequency  at  several  X/T  locations.  As  pointed  out  by  Devenport  and  Simpson 
(1990),  the  probability  density  function  of  the  switching  frequency  resembled  the 
log-normal  distribution  at  the  measuring  points  showing  a  bimodal  velocity  his¬ 
togram. 

Figure  55  shows  the  variation  of  the  average  switching  time  (T,u;)  in  the  bi¬ 
modal  zone.  The  average  switching  time  is  defined  as  the  average  period  for  the 

t 

flow  to  switch  back  to  the  previous  flow  mode.  For  Res  =  330  case,  the  average 
switching  times  {T^^)  at  X/T  =  -0.150  and  XjT  —  -0.175  were  4.72  and  4.43 
seconds,  respectively,  which  were  shorter  than  those  at  other  locations.  This  implies 
that  switching  from  one  mode  to  another  occurred  more  often  in  this  region.  As  the 
wing  was  approached,  the  mode  change  took  place  slowly  (T,u,  w  9(sec.))  especially 
at  low  Y  locations. 

For  the  Re$  =  1100  case,  the  switching  occurred  more  frequently  at  low  Y 
location  for  X/T  =  —0.125,  and  at  high  Y  location  for  X/T  =  —0.150. 

Figure  56  shows  the  variations  of  the  average  switching  frequency  {f^^)  in  the 
bimodal  zone  for  each  flow  condition.  The  switching  frequency  is  defined  as  the 
average  of  the  inverse  of  the  period  for  the  flow  to  switch  back  to  the  previous 
mode.  For  the  Reg  =  330  case,  the  values  of  the  switching  frequency  were  widely 
spread  across  the  bimodal  zone  ranging  from  0.20  Hz  to  0.80  Hz.  The  switching 
occurred  in  a  wide  Y  range  at  X  location  of  XjT  =  —0.150. 

For  the  Res  =  1100  case,  the  switching  frequency  varied  from  0.85  Hz  to  1.72 
Hz.  The  variations  of  across  the  bimodal  zone  at  each  X  location  explain 
the  intensity  of  the  switching  activity  as  follows.  At  XjT  =  -0.175,  which  was 
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the  most  upstream  location  showing  a  bimodal  histogram,  only  one  location  was 
observed  to  be  bimodal,  and  had  the  lower  switching  frequency  than  the  other 
points  at  different  X  locations.  The  variations  at  X/T  =  —0.15  show  that  the 
switching  appeared  to  occur  more  often  at  higher  Y  location  and  became  slower  at 
lower  Y  location.  On  the  other  hand,  this  trend  became  reversed  XfT  =  —0.125 
where  the  most  intense  large-scale  fluctuations  appeared.  The  averaged  switching 
frequency  at  X/T  =  -0.125  was  estimated  as  approximately  1.23  Hz.  As  the  wing 
was  approached  to  X/T  =  —0.10,  the  switching  activity  became  intense,  and  the 
values  of  the  switching  frequency  at  each  Y  location  were  approximately  0.89  Hz. 

Figure  57  illustrates  the  results  showing  the  fraction  of  the  backflow  mode  (74) 
during  each  point  measurement.  For  Re$  =  330  case,  the  back  flow  mode  was 
prevailing  for  most  of  the  time  near  the  wing  {X/T  =  —0.10  and  X/T  =  —0.125). 
However,  at  X/T  =  -0.15,  the  chances  were  almost  equal  for  each  flow  mode  to  take 
place  across  all  the  Y  locations.  For  the  flow  far  upstream,  the  zero  mode  flow  began 
to  last  longer  especially  near  the  wall.  The  same  trend  was  found  for  Ree  =  1100. 
However,  the  location  where  each  mode  took  place  at  approximately  equal  rate  was 
at  X/T  =  -0.125  near  the  wall,  which  means  that  the  frequent  switching  occurs 
near  the  wing  for  the  high  Reynolds  number  case.  Figure  58  shows  the  contour  plots 
of  the  backflow  mode  (76)  in  the  bimodal  zone.  The  dotted  lines  were  constructed 
based  on  the  assumption  that  74  does  not  vary  much  with  height  apart  from  the 
wall  (Devenport  and  Simpson,  1990a). 

The  switching  frequency  analysis  was  performed  on  an  IBM  PC/AT  compatible 
for  every  data  point  having  the  bimodal  velocity  histograms.  The  listing  of  the 
programs  developed  for  the  switching  frequency  analysis  is  tabulated  in  Table  7, 
and  a  diskette  containing  those  files  is  attached  on  the  back,  can  be  expressed 
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6.5  Results  of  PIDV 


The  three-dimensional  instantaneous  velocity  vector  field  was  obtained  in  the 
region  near  the  leading  edge  of  the  wing  from  the  high  speed  video  images. 

Figure  59  shows  a  sequence  of  the  bilevel  images  that  are  preprocessed  with 
the  techniques  including  the  backgrotmd  subtraction,  the  8-point  local  averaging 
and  the  thresholding.  The  images  in  the  left-hand  side  of  this  figure  represent 
side  views,  and  the  ones  in  the  right-handed  side  are  combined  images  for  the 
bottom  and  upstream  views.  The  combined  images  were  produced  by  the  mirror 
arrangement  described  in  Figure  15. 

Figure  60  shows  the  graphical  reconstruction  of  the  three-dimensional  instan¬ 
taneous  velocity  vector  field  in  the  nose  region  of  the  wing.  The  three-dimensional 
reconstruction  of  the  velocity  field  was  obtained  by  the  perspective  transformation 
from  the  three-dimensional  points  in  the  physical  plane  into  the  points  in  a  two- 
dimensional  plane.  In  addition,  animation  movie  wis  made  by  shooting  frame  by 
frame  using  Super  8mm  moviecorder.  The  animation  movie  shows  the  instantaneous 
velocity  field  information  and  the  trajectory  of  the  individual  bubble  elements. 

The  color  coding  was  used  to  represent  the  variations  in  Z-coordinate.  Since 
the  number  of  color  that  can  be  used  at  once  in  the  VGA  mode  of  the  maximum 
resolution  (640  x  480)  is  restricted  to  sixteen,  the  color  coding  looks  coarse. 

Figure  61  shows  the  three  instantaneous  velocity  components  of  the  individual 
hydrogen  bubble  elements.  This  was  obtained  by  the  3-D  PIDV  following  the 
individual  bubble  elements.  The  top  figure  represents  the  time  trajectory  of  the 
bubble  element,  and  the  bottom  one  contains  three  components  of  the  velocity  at 
the  corresponding  location.  Here,  the  hydrogen  bubble  elements  generated  from  the 
first  wire  are  plotted. 
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7.  DISCUSSIONS  OF  THE  RESULTS 


7.1  Description  of  the  Junction  Flow 

Flow  visualization  using  the  hydrogen  bubble  technique  revealed  details  of  the 
time-dependent  flow  structures  in  the  stagnation  region  of  the  wing,  which  was 
not  possible  by  the  surface  oil  flow  visualization  technique  or  single  point  LDV 
measurements  that  have  been  used  in  most  of  the  previous  studies. 

In  Figure  26,  at  the  instant  of  the  generation  of  the  inverted  time-line,  three 
distinct  re^ons  marked  by  “ui”,  and  “us”  were  observed  near  the  wall.  Each 
region  was  approximately  centered  at  JC/T  of  -0.4,  -0.25  and  -0.09,  respectively.  In 
the  outer  layer,  some  disturbance  as  seen  at  location  “si”  approached  the  wing. 
In  region  “V3”,  a  primary  junction  vortex  was  observed  while  a  weak  vortex  with 
the  same  rotation  direction  was  present  in  region  “V2”.  Occasionally,  a  vortex-like 
movement  was  observed  between  “U2”  and  “ua”,  which  will  be  explained  later  in 
this  section.  Near  the  bubble  generating  wire,  the  flow  seemed  to  lift  up  in  region 
“ui”.  The  cores  of  individual  vortices  are  indicated  by  “Ci”and  “C2”. 

In  Figure  26(b),  the  inverted  time-line  had  a  disturbed  region  as  seen  at  lo¬ 
cation  “d”.  In  Figure  26(b)  -  (h),  the  fluids  in  region  “ui”,  which  were  lifted  up, 
were  observed  to  leap-frog  over  the  larger  vortex  structure.  Here,  the  fluid  marked 
by  “d”  experienced  a  rotational  motion  emd  merged  into  the  primary  junction  vor¬ 
tex  from  the  outer  layer.  In  Figure  26(b)  -  (f),  the  two  vortex  structures  “Ci” 
and  “C2”‘were  apparent  although  some  disturbances  from  upstream  rolled  up  over 
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these  structures.  In  Figure  26(g),  large-scale  vortex-like  motions  indicated  by  “52” 
and  “53”  were  built  upstream  due  to  the  lifted-up  fluids  and  leap-frogged  over  two 
jimction  vortices,  C\  and  €2-  During  leap-frogging  of  the  small  vortices,  the  vortex 
marked  by  “C2”  lost  its  strength,  finally  merging  into  the  primary  jimction  vortex 
“Cl”,  as  shown  in  Figure  26(g)  -  (1). 

In  Figure  26(i)  -  (1),  the  flow  very  near  the  wall  was  divided  into  two  regions, 
“ui”  and  “ua”.  Region  “ui”  contained  low  momentum  fluids.  In  Figure  26(i),  some 
backflow  existed  in  region  “ui”.  In  Figure  26(j)  -  (1),  the  primary  vortex  moved 
downstream.  Later,  another  inverted  time-line  was  generated  in  Figure  26(1).  A 
fluid  region  in  the  outer  layer,  which  was  clearly  discernible  from  others,  is  marked 
by  “da”  In  Figure  26(m)  -  (u).  This  sequence  of  photographs  showed  that  a  vortical 
structure  existed  in  the  outer  layer,  which  did  not  merge  into  the  primary  junction 
vortex  upstream  of  the  wing.  However,  the  height  of  “da”  above  the  wall  reduced  as 
the  structure  approached  the  wing,  which  implies  that  it  might  be  affected  by  the 
primary  junction  vortex.  In  Figure  26(v),  it  was  not  possible  to  trace  “da”  since  it 
wrapped  around  the  wing. 

In  Figure  26(m)  -  (t),  the  lift-up  of  the  flow  was  again  observed  and  marked  by 

The  disturbances  created  in  the  outer  part  of  the  approaching  boundary  layer 
as  seen  at  location  “54”  in  Figure  26(o)  -  (r)  together  with  the  lifted-up  flow 
from  the- wall  layer  were  observed  to  merge  into  the  primary  junction  vortex.  Here, 
the  leap-frogging  of  fluid  regions  over  the  larger  vortical  structure  again  appeared. 

From  Figure  26(m)  on,  the  region  marked  by  “va”  appeared,  but  there  was  no 
distinct  vortex  structure  in  this  region.  Instead,  a  concentration  of  fluids  with  very 
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low  momenttun  seemed  to  exist  around  the  region  where  some  of  the  fluids 
merged  into  “ua”  and  some  followed  the  local  direction  of  the  lifted-up  fluid  flow. 

The  three  distinct  regions,  “ui”,  “^2”,  and  “U3”,  observed  from  the  flow  se¬ 
quence  observed  in  Figure  26,  correspond  to  the  lift-up  region,  the  secondary  vortex 
region,  and  primary  vortex  region,  respectively.  The  existence  of  these  regions  is 
consistent  with  the  LDV  measurements  as  in  Figure  43,  in  which  there  exist  three 
regions  with  higher  turbulent  stresses. 

The  sequence  of  the  photographs  (Figure  26)  further  shows  that  the  secondary 
vortex  region  merges  into  the  primary  vortex.  The  primary  vortex,  when  it  gained 
strength  to  some  extent,  moves  downstream  while  the  high  momentum  fluid  en¬ 
ter  this  region  from  upstream  just  like  the  “sweep-inrush”  phase  in  a  turbulent 
boundary  layer.  This  leads  to  the  stretching  of  the  primary  vortex. 

Flow  near  the  saddle  point  formed  between  the  primary  and  the  secondary 
vortex  was  captured  in  a  sequence  of  photographs  in  Figure  27.  In  this  sequence, 
fluid  region  (“T”)  followed  the  secondary  vortex  motion  marked  by  “S”  and  moved 
down  to  the  near  wall  region.  Later,  some  of  that  fluid  changed  the  flow  direction 
to  follow  the  primary  vortex  motion  (“P”).  This  shows  that  the  interaction  between 
two  active  vortices  created  a  tertiary  vortex  between  them  (see  Figure  27(j)).  This 
may  explmn  the  observation  of  the  multiple  vortex  system  existing  in  this  region 
by  the  previous  researchers. 

In  Figure  28,  the  leap-frogging  of  the  smaller-scale  structures  over  a  large  struc¬ 
ture  was  again  observed.  In  particular,  Figure  28(a)  -  (b)  illustrate  that  smaller- 
scale  eddies,  as  seen  at  locations  “ci”,  “02”,  “61”,  “62”)  and  “63”,  rolled  up  around  a 
primary  junction  structure  marked  by  region  “M”  to  merge  into  a  larger  structure. 
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In  Figure  28(c),  there  was  a  region  marked  by  “J”  where  the  distance  between  the 
time-lines  was  longer  than  the  others.  This  indicates  that  the  strong  activity  of 
the  junction  vortex  induced  the  flow  acceleration  in  this  region.  The  approximate 
location  of  the  core  of  the  primary  junction  vortex  is  marked  by  “C7”. 

Figure  28(d)  shows  an  instant  when  a  axial  component  of  vortex  filaments 
existed  at  location  “T”  and  started  to  rise  due  to  the  presence  of  the  junction  vortex. 
This  was  again  observed  in  Figure  28(e)  and  Figure  28(f).  The  region  marked  by 

in  Figure  28(e)  shows  that  the  flow  in  the  outer  region  was  accelerated  toward 
the  junction  vortex. 

In  Figure  28(c)  -  (g),  the  vortex  core  as  seen  at  location  “C”  was  observed  to 
change  its  location  upstream  of  the  wing  with  its  strength  varying  occasionally,  The 
strength  of  the  junction  vortex  may  be  estimated  by  the  intensity  of  the  white  area 
in  the  picture  formed  by  the  agglomeration  of  the  hydrogen  bubbles.  The  region 
consisting  of  low  momentum  fluids  was  observed  upstream  of  the  junction  vortex 
although  this  was  hardly  discernible  in  the  photographs.  In  most  of  the  pictures, 
the  low  momentum  fluid  region  existed  between  the  wire  location  and  the  primary 
junction  vortex. 

Figure  28(h)  shows  an  instant  when  the  time-lines  proceeded  to  the  leading 
edge  of  the  wing  without  much  interference,  which  may  indicate  the  junction  vortex 
was  stretched  around  the  nose  of  the  wing  as  seen  at  the  region  “5/”.  Indeed,  the 
flow  near  the  wall  was  accelerated,  and  the  vortex  core  indicated  by  the  white  area 
was  not  observed  due  to  the  effect  of  the  stretching.  It  should  be  notable  that 
the  instability  observed  along  the  leading  edge  is  possibly  a  result  of  the  amplified 
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unsteadiness  in  the  vicinity  of  the  stagnation  point  due  to  the  high  free-stream 
turbulence  level. 

Figure  28(i)  again  shows  a  axial  component  of  vortex  filaments  formed  upstream 
of  the  jvmction  vortex  which  was  observed  in  re^on  “il”. 

The  location  marked  by  '‘52”  in  Figure  28(j)  shows  a  large  structure  in  the 
outer  layer.  The  time-lines  surrounding  this  structure  indicate  that  the  convective 
speed  of  the  large  structure  in  the  outer  layer  is  slower  than  the  velocity  of  the 
external  flow.  However,  the  faster  speed  near  the  wall  region  indicates  that  the 
vortex  was  being  stretched. 

Figure  28(k)  shows  that  the  multiple  vortex-like  disturbances  around  a  large 
vortical  structure  near  the  junction  contributed  to  the  formation  of  the  primary 
vortex  structure.  The  location  marked  by  “fci”  shows  the  engulfment  from  the 
outer  energetic  flow  into  the  region  of  low  momentum  fluid  flow  to  make  a  stronger 
vortical  structure. 

As  observed  from  flow  visualization,  a  repetitive  flow  behavior  was  found  in 
this  region,  and  may  be  considered  to  be  made  up  of  several  characteristic  events. 
Yet,  on  some  occasions,  some  of  the  events  could  not  be  observed.  Those  events 
observed  from  flow  visualization  can  be  summarized  as  follows: 

i  )  Primary  and  secondary  vortices  are  formed  upstream  of  the  wing. 
it  )  The  backflow  leads  to  the  lift-up  of  the  fluid  in  the  region  of  low  momentum 
flow.  The  backflow  is  produced  in  the  mean  sense  by  the  formation  of  the 
secondary  vortex  upstream  of  the  primary  vortex.  This  lift-up  of  fluid  generates 
the  disturbances  in  the  flow. 
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t«  )  Multiple  vortices  exist  in  the  form  of  disturbances  or  small  structures,  up¬ 
stream  of  the  primary  junction  vortex.  In  some  instances,  the  axial  vortex 
filaments  are  observed. 

iv  )  Smaller-scale  structures  leap-frog  over  a  primary  junction  vortex  to  form  a 
stronger  vortical  structure.  The  secondary  vortex  also  merges  into  the  primary 
vortex. 

V  )  When  the  primary  junction  vortex  gzuns  sufficient  strength,  it  seems  to  stretch 
at  a  faster  rate  around  the  nose  of  the  wing,  and  this  was  observed  repetitively. 

vi  )  Very  unstable  flow  follows  after  faster  stretching  while  flow  is  stable  during 
the  stretching  due  to  the  activity  of  the  primary  junction  vortex. 

vii  )  After  a  while,  the  backflow  reaches  far  upstream  to  make  a  large,  but  weak 
junction  vortex  structure.  Here,  the  primary  and  secondary  vortex  are  formed 
again. 

The  above  sequence  of  events  is  described  in  Figure  62.  The  percentage  at  the  end 
of  each  comment  represents  the  time  proportion  of  that  event  in  the  average  period 
of  the  acyclic  flow. 

This  acyclic  process  helps  expleiin  the  existence  of  the  two  flow  modes  suggested 
by  Devenport  and  Simpson  (1987).  When  the  secondary  vortex  exists,  the  backflow 
reaches  far  upstream,  which  appears  on  the  average  as  the  backflow  mode.  The  zero 
flow  mode  can  be  explained  eis  the  merging  of  the  secondary  vortex  into  the  primary 
one,  and  the  stretching  of  the  primary  vortex  accompanying  the  entrainment  of  the 
high  momentum  fluid  into  this  region.  The  acyclic  behavior  of  the  junction  flow 
seems  to  be  driven  by  the  repetitive  increase  in  the  rate  of  stretching  of  the  vortex 
system.  ’ 
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Some  statistics  on  the  flow  visualization  restdts 


Some  statistical  information  was  extracted  from  the  flow  visualization  by  in¬ 
specting  the  specific  points  in  the  flow  region  or  specific  events  in  the  flow.  Repeti¬ 
tive  measurements  were  possible  by  watching  the  recorded  tapes  as  much  as  needed. 
Ten  minutes  of  the  flow  visualization  results  was  used  to  get  quantitative  data  as 
shown  in  Table  9. 

The  switching  frequency  of  two  flow  modes  was  studied  apart  from  the  wall 
at  X/T  —  -0.15  by  counting  time  for  each  switching  to  occur.  The  switching 
was  decided  by  visual  inspection.  The  fraction  of  a  backflow  mode  (7^)  at  the 
same  location  was  obtained  statistically  by  counting  the  total  duration  time  in  the 
backflow  mode.  From  ten  minutes  of  the  flow  visualization  videos,  the  average 
switching  frequency  (/,„,)  was  approximately  0.18(Hz),  and  74  was  approximately 
0.78.  Compjiring  to  the  results  from  SFA  in  Figure  57,  the  time  fraction  of  the 
backflow  mode  (74)  is  considered  to  be  overestimated.  The  corresponding  Strouhal 
number  (5tj  =  f sruT/Unf)  was  approximately  0.28. 

The  average  lapse  time  between  the  formation  of  successive  structures  was 
estimated  to  be  around  1.26±0.39(sec).  This  number  represents  the  average  period 
of  the  acyclic  flow  structure.  The  average  time  proportion  for  a  primary  vortex  to 
last  before  the  separation  vortex  forms  was  about  20%  of  the  total  period,  which 
is  about  0.25  ±  0.08(sec).  The  average  time  for  the  leap-frogging  process  to  occur 
was  approximately  0.55  ±  0.16(sec),  which  corresponds  to  44%  of  the  total  period. 
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7.2  Estimation  of  friction  velocity  (ur) 


The  friction  velocity  (ur)  was  estimated  and  compared  by  two  methods  :  the 
“cross  plot  of  log-region”  procedure  used  by  Rvmstadler  et  al.  (1963)  and  the  “wall- 
profile  method”  by  Schraub  and  Kline  (1965).  The  procedures  and  results  of  the 
two  methods  are  described  in  Figure  32  -  Figure  37.  It  should  be  noted  that  since 
the  procedure  by  the  “cross  plot  of  log-region”  assumes  the  universal  law  of  the 
wall,  the  resulting  depends  on  which  part  of  data  is  forced  to  fit.  The  current 
study  adopted  the  universal  logarithmic  relation  given  by  Clauser  (1954)  : 

u"*"  =  2.44  In  +  4.9  (7.1) 

Therefore,  this  method  becomes  unreliable  if  there  are  regions  near  the  wall  where 
the  law  of'the  wall  is  not  universal. 

The  wall-profile  method  assumes  the  existence  of  a  linear  region  where  u'^  —  y'*’ 
holds.  Thus,  this  method  to  estimate  a  value  of  wall  shear  can  be  applied  if  there 
are  enough  data  points  in  the  viscous  sublayer  region  for  fitting. 

Figure  32  shows  the  repetitive  procedure  to  obtain  the  local  skin  friction  coeffi¬ 
cient  for  the  high  Reynolds  number  case  by  forcing  the  data  points  in  the  logarithmic 
region  to  the  adopted  straight  line.  The  local  skin  friction  coefficient  of  0.0053  was 
obtained  by  this  method.  With  the  friction  velocity  calculated  using  this  estimated 
skin  friction  coefficient,  the  law  of  the  wall  plot  was  made  in  Figure  33.  Some 
deviation  of  the  data  points  from  the  relation  («■*■  =  y'^)  was  observed. 

However,  the  “wall  profile”  method  used  in  Figure  34  suggested  that  another 
value  of  local  skin  friction  coefficient  (C/  =  0.0078)  works  well  in  the  viscous  sub¬ 
layer  region,  and  the  universal  relation  in  the  logarithmic  region  might  have  the 
different  coefficients  (it"*"  =  2.44 In y"*"  -f- 1.9). 
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Figure  35  -  37  shows  the  same  procedure  for  the  low  Reynolds  number  case. 
The  local  skin  friction  coefficients  of  0.0070  and  0.0077  were  obtained,  respectively. 
In  this  case,  the  logarithmic  relation  did  not  deviate  much  from  the  Clauser  coef¬ 
ficients,  but  some  differences  were  again  observed.  In  either  test  condition,  some 
deviations  in  the  viscous  sublayer  occurred  when  the  “cross  plot  of  the  log-region” 
procedure  was  used.  On  the  other  hand,  the  “w£dl  profile”  method  created  devia¬ 
tions  in  the  log  region. 

In  general,  the  resulting  wall  shear  estimated  using  each  method  does  not 
coincide  very  well  as  pointed  out  by  Schraub  et  al.  Table  1  presents  a  summary  of  the 
estimated  friction  velocities.  Differences  up  to  18%  existed  between  u^’s  determined 
by  the  two  methods.  The  uncertainty  of  determining  the  friction  velocity  for  each 
method  separately  is  estimated  to  be  ±7%  as  shown  in  Table  2. 

7.3  Mean  Velocity  Profiles 

The  mean  location  of  the  core  of  the  primary  junction  vortex  was  estimated 
from  the  mean  velocity  profile  plot  by  locating  the  intersection  between  the  maxi¬ 
mum  gradient  of  the  mean  velocity  and  the  zero  mean  velocity.  The  approximate  lo¬ 
cations  were  at  X/T  =  -0.125  and  Y/T  =  0.032  for  Reg  =  1100,  and  X/T  =  -0.150 
and  Y/T  =  0.04  for  Ree  =  330.  Comparioon  with  Devenport’s  data  {X/T  =  —0.20 
and  Y/T  =  0.05)  indicates  that  the  locations  of  the  vortex  cores  did  not  show 
good  correlation  when  scaled  with  T,  but  the  physical  .iV-location  approached  the 
leading  edge  of  the  wing  and  the  physical  Y -location  weis  getting  lowered  with  the 
increasing  Reg.  This  implies  that  the  mean  location  data  for  the  vortex  cores  may 
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be  correlated  when  different  nondimensionalizing  parameters  aie  used  for  scaling. 
For  instance,  the  use  of  6  gives  X/6  =  —0.530  and  Y/6  =  0.141  for  Ree  =  330,  and 
X/6  =  —0.442  and  Y/6  =  0.113  for  Ree  =  1100.  Devenport’s  data  can  be  scaled  to 
X/6  =  —0.4  and  Y/6  =  0.10  (i2e#  =  6700).  These  nondimensional  mean  locations 
of  the  vortex  cores  for  each  case  show  some  Reynolds  number  dependency  that  the 
vortex  is  getting  closer  to  the  junction  with  the  increasing  Reynolds  number. 

It  is  considered  that  the  Reynolds  number  did  not  change  the  mean  flow  be¬ 
havior  qualitatively.  However,  some  of  the  plots  in  Figure  38  and  Figure  39  show 
that  there  existed  very  similar  profiles  among  three  momentum  thickness  Reynolds 
number  cases.  As  with  the  case  of  the  mean  vortex  core  location,  the  use  of  ^  as  a 
scaling  parameter  gives  a  correlation. 

7.4  Low  Frequency  Unsteadiness 

Figure  63  shows  the  typical  LDV  signals  observed  in  the  region  where  the 
backflow  was  found  to  be  very  dominant.  As  pointed  out  by  Devenport  and  Simpson 
(1988a,  1990),  it  seemed  like  there  existed  two  modes  of  flow  (a  backflow  mode  and 
a  zero  flow  mode),  and  these  two  modes  of  flow  were  switching  back  and  forth  in  the 
vicinity  of  the  vortex  core.  In  addition,  the  evidence  of  flow  visualization  reveals  that 
the  flow  in  this  region  has  an  “acyclic”  behavior.  The  vortex  system  was  observed 
to  move  around  the  nose  region  of  the  wing  with  low  frequency  fluctuations. 

In  this  viewpoint,  the  flow  velocity  may  be  considered  to  consist  of  three  parts, 
which  are  a  mean  velocity  term,  a  large-scale  fluctuation  term,  and  a  small-scale 
fluctuation  term.  The  mean  velocity  term  and  small-scale  fluctuations  are  regarded 
as  the  Reynolds  decomposition  terms.  On  the  other  hand,  the  third  term,  i.e.,  the 
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large-scale  fluctuation  term,  is  superimposed  on  this  turbulent  flow  at  low  frequency 
due  to  the  presence  of  the  wing-body  junction.  Three  component  decomposition  of 
the  flow  velocity  can  be  written  as  follows: 

U=^U  +  u',  +  u'i.  (7.2) 

These  large-scale  fluctuations  (u{)  are  considered  to  be  responsible  for  the  large- 
scale  time-dependent  junction  flow  activities  such  as  the  bimodal  flow  structure 
found  by  Devenport  and  Simpson  (1986, 1987,  1988a,  1988b,  1990). 

The  general  tendency  in  u'^  plots  reveals  that  the  turbulent  activity  is  more 
intense  at  high  Y  locations  in  the  upstream  flow  than  at  low  Y  location,  and,  as 
the  flow  moves  downstream,  the  intense  turbulent  region  moves  down  near  the 
wall.  This  trend  was  also  found  in  the  switching  frequency  analysis  discussed  in  the 
previous  chapter.  The  switching  between  the  two  flow  modes  occurred  more  often 
at  high  Y  location  upstream  and  at  low  Y  near  the  wing. 

It  should  be  noted  that  the  similar  behavior  of  u'^  variations  in  three  different 
Reynolds  number  cases  indicates  the  possibility  of  the  existence  of  other  scaling 
parameters.  Here,  the  boundary  layer  thickness  {6)  seems  to  give  better  scaling 
between  the  three  cases. 

Figure  43  shows  the  contour  plots  of  u'^  for  Re$  =  330  and  Ree  =  1100, 
respectively.  In  Figure  43(a),  the  most  intense  turbulent  region  was  found  very 
near  the  wall  approximately  at  X/T  =  —0.175  and  another  distinct  region  was 
observed  approximately  at  X/T  =  —0.3.  Upstream  of  the  main  junction  vortex 
structure,  which  can  be  identified  by  the  intense  turbulent  region,  the  region  of  fluid 
consisting  of  low  turbulent  intensity  weis  observed.  This  region  is  considered  eis  the 
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one  conesponding  to  region  of  the  low  momentum  fluids,  which  was  discussed 
in  the  previous  section.  For  the  high  Re$  case  as  shown  in  Figure  43(b),  the  vortex 
structure  which  contained  a  large  amount  of  turbulent  energy  was  confined  near 
the  leading  edge  of  the  wing  (X/T  =  —0.125).  Here,  the  small  region  of  the  low 
turbulence  intensity  was  again  found  upstream  of  the  primary  structure. 

The  results  obtained  from  the  analysis  of  the  switching  frequency  in  the  bimodal 
region  (Figure  56)  show  that  the  switdiing  frequency  ranges  from  0.2Hz  to  0.8Hz  for 
Res  =  330,  and  from  0.85Hz  to  1.72Hz  for  Res  =  1100.  Schwind  (1962)  observed 
two  types  of  unsteady  oscillation  according  to  the  flow  velocity  when  he  studied  a 
flow  around  a  60®  included-angle  wedge.  For  the  flows  of  regime  4  type  (Schwind ’s 
definition),  a  frequency  of  oscillation  of  0.8Hz  was  measured,  and  a  frequency  of 
1.4Hz  for  the  flows  of  regime  5  type.  The  switching  frequency  analyzed  in  this 
study  and  the  frequency  of  oscillation  measured  by  Schwind  represent  different 
characteristics  of  flow.  Yet,  it  could  seem  likely  that  the  low  speed  case  {Ree  = 
330)  may  be  considered  as  Schwind ’s  regime  4  or  5,  and  the  higher  speed  case 
{Res  =  1100)  as  Schwind’s  regime  5  or  above. 

7.5  Discussion  of  higher  order  turbulence  results 

The  skewness  and  flatness  factors  for  the  U-component  velocity  fluctuations 
were  obtained  at  all  the  merisurement  locations  from  the  present  LDV  data.  Fig¬ 
ure  45  shows  the  flatness  factor  distributions  for  the  Res  =  1100  case.  At  X/T  = 
—0.86,  the  velocity  distribution  is  still  close  to  that  of  the  Gaussian  (Fu  =  3)  in 
most  of  the  outer  layer.  As  the  wing  is  approached  but  still  before  separation,  the 
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flatness  factors  become  deviated  from  the  value  of  the  Gaussian  distribution  near 
the  wall  and  in  the  outer  part  of  the  shear  layer,  while  most  of  measurement  point 
in  the  shear  layer  has  the  velocity  distribution  still  similar  to  the  normal  distribu¬ 
tion.  In  the  vortical  flow  region,  the  flatness  factors  are  far  from  being  the  normal 
distribution. 

Figure  46  shows  the  skewness  factor  distribution  plots.  As  pointed  out  by  Simp¬ 
son,  Chew,  and  Shivaprasad  (1981),  downstream  of  separation  the  skewness  factor 
for  V,  becomes  negative  values  near  the  wall,  while  large  positive  skewness  factors  are 
found  upstream  of  separation.  This  is  largely  due  to  the  reduced  large-amplitude 
fluctuations  in  the  backflow  region.  The  skewness  factors  show  the  tendency  to 
change  signs  at  the  location  of  maximum  Reynolds  stresses  or  large  fluctuations. 
The  intense  exchange  of  momentum  in  this  region  reduces  the  very  large  or  very 
small  fluctuations  of  u  component  velocity.  In  the  outer  part  of  the  shear  layer,  the 
skewness  factor  is  negative  for  almost  all  the  measurement  points. 

For  the  Reg  =  330  case  as  shown  in  Figure  47  and  Figure  48,  the  trend  of  both 
factors  shows  similar  behavior  across  the  shear  layer  at  all  the  measurement  points. 

7.6  Dimensional  Analysis  on 

A  dimensional  analysis  on  the  switching  frequency  was  made  to  obtain  the  di¬ 
mensionless  parameters  associated  with  the  low  frequency  fluctuations.  The  average 
switching  frequency,  can  be  expressed  as  the  following  functional  form  with 
the  other  related  parameters  for  a  given  cylindrical  wing  shape; 

Jsw  =  9{Uref,  T,a,X,  Y,  6,  d,  H,  p,  p)  (7.3) 
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where  is  the  average  switching  frequency,  Uref  is  the  reference  velocity  measured 
at  X/T  =  —5.312,  T  is  the  maximum  thickness  of  the  wing,  a  is  the  length  of  the 
major  axis  of  the  elliptical  nose,  X  and  Y  are  the  locations  upstream  of  the  wing,  6  is 
the  thickness  of  the  approach  boundary  layer  at  X/T  =  -2.15,  B  is  the  momentum 
thickness  at  the  same  location,  H  is  the  shape  factor,  and  fi  and  p  are  the  fluid 
viscosity  and  density,  respectively. 

The  dimensional  analysis  on  Equation  (7.3)  gives 

where  Sir  =  77^  and  Rer  = 

Uref 

Here,  5t j-  is  the  Strouhal  number  associated  with  the  low  frequency  fluctua¬ 
tions,  and  Rex  is  the  Reynolds  number  based  on  the  maximum  thickness  of  the 
wing.  Since  only  one  shape  of  the  model  was  investigated  in  this  study,  the  depen¬ 
dence  on  the  slenderness  ratio  can  not  be  checked  and  this  dimensionless  parameter 
was  left  out  from  the  functional  form.  The  shape  factor,  H,  did  not  show  much  vari¬ 
ation  as  seen  in  Table  1,  and  was  assumed  not  to  affect  the  dependent  dimensionless 
variables.  Thus,  the  remaining  terms  can  be  written  as 

^  X  Y  6, 

Six  g\RcT ,  y  ’  ) 

or  (7.5) 
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This  means  that  the  switching  frequency  is  represented  by  the  Strouhal  number 
(5tT)»  aiid  this  Strouhal  number  is  a  function  of  the  Reynolds  number  based  on  T, 
the  measurement  location  upstream  of  the  wing  (XfT,YfT),  and  the  characteristics 
of  the  approach  boundary  layer  (6/T  or  d/T). 

Figure  64  shows  the  variations  of  the  Strouhal  numbers  at  various  locations 
upstream  of  the  wing  for  each  flow  condition.  For  the  approximate  location  of 
the  vortex  core  {X/T  =  -0.125)  the  Strouhal  number  showed  near  constant  value 
of  0.S4  for  Reg  =  330,  and  0.40  for  Reg  =  1100.  For  the  Reg  =  330  case,  the 
Strouhal  number  agreed  fairly  well  with  the  result  obtained  from  flow  visualization 
{Sir  =  0.28  from  Table  9).  For  the  Reg  =  1100  case,  the  Strouhal  numbers  across 
the  whole  bimodal  zone  did  not  vary  much,  while  the  Strouhal  numbers  showed 
very  much  variation  at  all  the  locations  for  Reg  =  330. 

These  indicate  that  at  the  specific  Reynolds  number  the  Strouhal  number  based 
on  the  maximum  thickness  of  the  wing  may  be  a  fairly  constant  with  Y  position 
near  the  mean  location  of  vortex  core.  These  results  are  in  good  agreement  with 
the  spectral  analyses  made  on  the  hot  wire  signals  measured  by  Baker  (1979).  His 
measurements  were  made  at  the  height  of  1mm  from  the  wind  tunnel  wall  and  at 
X/D  =  -0.82  which  were  outside  of  the  primary  vortex  region  but  still  in  the  sep¬ 
arated  region.  His  results  showed  that  Sip  =  0.26  at  low  speed  (  D/6*  =  33.8  ), 
and  Sip  —  0.40  at  high  speed  (  D/6*  =  34.2  ).  He  found  that  the  high  frequency 
oscillation  Stp  increases  as  UD/u  increases  {Stp  increasing  up  to  about  0.6  at 
higher  UD/u).  Devenport  and  Simpson  (1990a)  obtained  switching  frequencies  at 
X/T  =  -0.20  by  constructing  the  histograms  and  considering  bimodal  histograms 
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as  formed  by  two  added  Gaussian  lustograms.  At  Uref  =  21ml s,  they  got  the  aver¬ 
age  frequency  of  38  Hz  and  the  most  common  frequency  of  91  Hz,  which  correspond 
to  the  Strouhal  number  of  0.1  and  0.24,  respectively.  Olgmen  et  al  (1991)  made 
another  switching  frequency  analysis  of  the  pressure  fluctuation  data  obtained  near 
the  vortex  core  in  the  wind  timnel  at  Urtf  =  32.5m/s.  In  his  measurements,  the 
half-size  model  of  the  same  geometry  was  used  and  the  different  reduction  method 
was  applied  to  obtain  the  switching  frequency.  He  obtained  the  switching  frequency 
of  68  (Hz),  which  corresponds  to  the  Strouhal  number  of  0.15. 

Several  dimensionless  frequency  groups  were  tried  using  different  parameters 
for  scaling  between  several  data  sets  :  i)  Unf/T,  it)  Uref/^,  and  Hi)  {dU lds)noaf 
As  will  be  described  in  the  next  section,  the  value  in  tit)  represents  the  rate  of 
vortex  stretching  and  was  obtained  by  an  inviscid  panel  method  (Fleming,  1991). 
Comparison  of  the  results  obtained  for  several  sets  of  data  are  provided  in  Table  10. 
Among  those  groups,  6  gives  least  variation  of  the  results  between  data  sets. 

7.7  Vortex  Stretching 

Combined  with  the  descriptive  model  for  the  acyclic  process  observed  from 
flow  visualization,  the  switching  frequency  between  two  flow  modes  discussed  in  the 
previous  section  heis  a  direct  relation  with  the  vortex  stretching.  As  discussed  in 
the  section  7.1,  the  vortex  stretching  occurs  almost  every  cycle,  and  this  creates 
the  inrush  of  high  momentum  fluid  into  the  separated  flow  region  (zero  flow  mode), 
while  the  secondary  vortex  contributed  to  forming  the  backflow  mode. 

Therefore,  the  switching  frequency  is  approximately  equal  to  the  frequency  of 
the  acyclic  process  observed  in  the  junction  region.  This  implies  that  the  Strouhal 
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number  ^opted  in  the  cvnrent  study  is  directly  related  to  the  rate  of  the  vortex 
stretching. 

The  vorticity  transport  equation  can  be  written  as  follows: 

^  =  a.  ^V  +  (7.6) 

which  shows  that  the  rate  of  change  of  vorticity  is  equal  to  the  rate  of  deformation 
of  the  vortex  lines  (vortex  stretching)  plus  viscous  diffusion  of  the  vortex. 

Here,  the  vortex  stretching  term  can  be  rewritten  as 


12  •  VV  = 


(7.7) 


where  tf,  is  a  unit  vector  in  streamwise  direction.  For  the  same  geometry, 


dV  Uoo 
ds  ^  T  ' 


(7.8) 


which  indicates  that  the  rate  of  the  vortex  stretching  depends  largely  on  the  flow 
velocity  (boundary  layer  characteristics). 

Therefore,  the  Strouhal  number  based  on  model  thickness  can  be  related  to  the 
vortex  stretching  as  follows: 


/ 

dV\ 


(7.9) 


This  implies  that  for  a  constant  Strouhal  number,  the  switching  frequency  is  linearly 
varying  with  the  rate  of  the  vortex  stretching.  It  is  noted  that,  for  the  same  amount 
of  vortex  strength,  the  rate  of  vortex  stretching  was  much  higher  in  the  wind  tunnel 
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measurements  (JRc#  =  4450)  than  in  the  water  tunnel  measurements.  They  are 
compared  as  follows: 

■ 

~  120  for  Ree  =  1100 

\~5s  /tuoter 

-^*  <-**'  -  595  for  Ree  =  330 

Considering  this  large  difference  in  the  stretching  rate  of  vorticity  between  air  and 
water,  the  results  obtained  from  the  switching  frequency  analysis  in  the  previous 
section  show  good  correlation  between  two  test  cases  {Sir  =  0.40  for  Re$  =  1100 
in  water  and  Sir  =  0.15  for  Re^  =  4450  in  air  by  Oilmen  ei  al). 

7.8  3-D  PIDV  Results 

The  sequence  of  the  bilevel  images  as  shown  in  Figure  59  illustrates  lots  of 
practical  difficulties  in  implementing  the  motion  analysis.  First  of  all,  the  gray  levels 
of  the  background  image  were  not  uniform  from  freime  to  frame.  This  fluctuation  of 
gray  levels  could  not  be  removed  because  it  was  caused  by  the  limitations  of  the  high 
speed  video  system  employed.  Therefore,  the  subtraction  of  the  same  background 
level  fium  all  the  frames  may  cause  the  problem  of  losing  information.  However, 
since  the  processing  speed  is  greatly  improved  by  employing  a  uniform  background 
subtraction,  some  routines  were  added  to  retrieve  the  lost  information.  In  order  to 
compensate  for  the  gray  level  fluctuations,  several  frames  (previous  10  frames)  were 
considered  to  determine  the  matching  between  the  bubble  elements.  Thus,  even 
though  we  lose  some  bubble  elements  at  some  instant,  they  may  be  retrieved  after 
several  frames. 
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Secondly,  the  bubble  overlapping  (see  Figure  24)  was  the  most  troublesome 
problem  in  processing  the  image  data.  Due  to  the  overlapping  of  the  bubble  el¬ 
ements,  the  information  of  at  least  one  element  was  lost  and  the  uncertainty  in 
determining  the  centroid  locations  was  increased  by  the  factor  of  two.  To  improve 
this  problem,  the  split-merging  detection  algorithm  (see  Figure  25)  was  developed. 
In  the  sequence  of  instantaneous  velocity  vector  field  as  shown  in  Figure  60,  some 
of  the  bubble  elements  were  retrieved  using  this  algorithm. 

As  seen  in  Figure  60,  most  of  the  bubble  elements  were  lost  near  the  wing  and 
in  the  vortex  region  due  to  the  fast  diffusion  of  the  hydrogen  bubbles.  This  is  the 
inherent  problem  which  cannot  be  avoided  when  the  hydrogen  bubbles  are  used  as 
tracer  particles. 

Figure  61  shows  the  trajectories  of  the  individual  bubble  elements  and  the 
corresponding  three-dimensional  instantaneous  velocity  components.  It  is  notable 
that  axial  acceleration  exists  at  t  =  0.35(sec.)  for  all  the  bubble  elements,  which 
may  be  explained  as  the  flow  event  when  the  vortex  is  stretched.  Among  the  three 
velocity  components,  Z-component  fluctuation  is  dominant  at  low  frequency. 

As  shown  in  Figure  60  and  Figure  61,  the  five-wire  hydrogen  bubble  tech¬ 
nique  presented  a  new  method  to  obtain  three-dimensional  velocity  components. 
However,  due  to  the  limitation  of  the  bubble  patch  size,  the  results  were  tilways 
subject  to  the  sizable  uncertainties.  Another  important  factors  that  affected  the 
results  were  fast  diffusion  of  hydrogen  bubbles  in  a  vortical  flow  and  distortion  of 
the  bubble  patches.  In  order  to  reduce  those  effects  in  a  highly  distorted  flow,  it 
is  recommended  to  use  another  tracking  particles,  which  are  small  in  size  and  not 
distorted  in  a  shear  layer  region. 
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8.  CONCLUSIONS 


Experimental  investigations  of  the  flow  structures  in  the  nose  region  of  a  wing- 
body  junction  were  made  in  the  water  tunnel  at  low  speeds. 

Several  flow  visualization  techniques  were  developed  using  the  hydrogen  bubble 
method.  An  inverted  time-line  pattern  allows  the  easier  tracking  of  the  inverted 
time-lines  while  the  whole  flow  fleld  is  visualized  by  a  sheet  of  hydrogen  bubbles. 
A  sequenced  multiple-wire  hydrogen  bubble  technique  was  developed  and  used  to 
visualize  the  three-dimensional  flow  field. 

LDV  measurements  were  performed  in  the  plane  of  symmetry  upstream  of  the 
wing.  The  boundary  layer  characteristics  obtained  at  XjT  =  —2.15  were  studied. 
The  mean  location  of  the  core  of  the  primary  vortex  is  more  confined  to  the  junction 
with  the  increasing  Reynolds  number.  Scaling  the  Y  coordinates  with  the  approach 
boundary  layer  thickness  gives  better  correlation  between  the  test  cases  of  different 
Reynolds  number  than  scaling  with  model  thickness. 

A  visual  study  reveals  that  there  exists  an  acyclic  flow  behavior  in  this  region, 
which  is  closely  related  to  the  stretching  of  the  horseshoe  vortex  system.  Three 
component  decomposition  of  the  flow  velocity  was  attempted  to  expleiin  the  large- 
scale,  low-frequency  fluctuations  present  in  this  region. 

The  bimodal  velocity  histograms  found  in  the  wind  tunnel  experiments  by  De- 
venport  and  Simpson  (1987,  1988a,  1988b,  1990)  were  also  obtained  iii  this  study. 
The  idea  of  switching  between  two  flow  modes  was  implemented  to  calculate  the 
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frequeiu^  of  switching  in  a  dimensionless  form.  The  Strouhal  number  at  the  ap¬ 
proximate  core  of  the  primary  vortex  represents  the  frequency  of  switching,  and  was 
obtained  as  follows; 

StT  =  0.34  for  Re$  =  330 

Str  =  0.40  for  Re$  =  1100. 

The  Strouhal  number  obtained  from  flow  visualization  results  {Six  =  0.28)  showed 
good  agreement  with  the  results  from  LDV  data. 

The  particle  image  displacement  velocimetry  (PIDV)  using  a  multiple  wire 
hydrogen  bubble  technique  was  developed  to  obtain  the  three-dimensional  instan¬ 
taneous  velocity  field  in  the  junction  region.  New  algorithms  to  implement  image 
processing  technique  and  motion  analysis  were  successfully  developed.  However, 
the  application  of  the  present  PIDV  technique  to  the  hydrogen  bubble  images  was 
improper  especially  in  a  highly  distorted  flow  region  due  to  the  fast  diffusion  of  hy¬ 
drogen  bubbles.  Therefore,  the  use  of  another  tracking  media  such  as  solid  particles 
is  recommended  in  a  shear  flow  application  of  the  present  PIDV  technique. 
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TABLES 
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I 


Table  1.  Test  conditions  at  X/T  =  -2.15 


Boundary  layer  characteristics 


case 

high  speed 

low  speed 

UrtAcm/s) 

44.75 

9.07 

S  (cm) 

2;73 

4.06 

6  (cm) 

0.254 

0.381 

H 

1.306 

1.430 

Cf^ 

0.0053 

0.007 

Ur^  (cm/s) 

2.155 

0.495 

0.0078 

0.0077 

Ur^  (cm/s) 

2.615 

0.518 

Ree 

1100 

330 

*  (measured  sA  X/T  —  —5.312) 

^  (estimated  by  the  “cross  plot  of  log-region”  procedure) 

*  (estimated  by  the  “wall  profile”  method) 


Table  2.  Uncertainties  of  measured  data 


Uncertainties 


Velocity  (LDV) 

U±3% 

Velocity  fluctuations 

u>^±4% 

FVinge  spacing 

df±  0.8% 

Probe  volume  location 

y  ±  2.5  X  10-*mm 

Friction  velocity 

Ur  ±  7% 

Hydrogen  bubble  probe  location 

X  ±  1mm 

Velocity  (PIDV) 

U  ±  7.5% 

Table  3.  Camera  calibration  coefficients 


Imager  I  (Side  view) 


68.5751 

0.6138 

-5.6227 

230.1952 

-0.0181 

68.7164 

-0.0913 

14.2271 

0,0010 

0.0009 

-0.0334 

1.0 

Imager  II  (Bottom  view) 


-0.7128 

38.5617 

0.4224 

112.9080 

-0.1299 

-0.0217 

37.9156 

93.6505 

-0.0143 

0.0005 

0.0022 

1.0 

Table  4.  Camera  characteristics 


Imager  resolution 

240  X  192  (lines) 

Field  of  view 

8.13(cm)  X  7.11(cm) 

Estimated  maximum  axial  velocity 

12.7(cm/s.) 

Frame  rate 

l/125(scc.) 

1  pixel  distance  in  axial  direction 

0.034(  cm) 

Estimated  movement 

of  centroid  location 

between  successive  frames 

0.102(cm) 

Estimated  pixel  distance 
of  centroid  location 

between  successive  frames 

3  (pixels) 

Table  5.  Test  log  for  LDV  measurements 


Profile  # 


14 


15 


16 


17 


Re$ 

XjT 

n 

B9 

330 

-0.86 

49870.0 

76290 

1100 

-0.86 

49870.0 

76290 

330 

-0.67 

49870.0 

76290 

330 

-0.46 

49870.0 

76290 

330 

-0.40 

49870.0 

76290 

330 

-0.35 

49870.0 

76290 

330 

-0.30 

49870.0 

76290 

330 

-0.25 

i  49870.0 

76290 

330 

-0.20 

49870.0 

76290 

330 

-0.15 

49845.2 

76290 

330 

1 

p 

o 

49826.9 

76290 

330 

-0.05 

49866.1 

76290 

330 

-0.025 

49870.0 

76290 

1100 

-0.025 

53063.7 

76290 

1100 

-0.05 

52853.2 

76290 

1100 

-0.10 

51497.4 

76290 

1100 

-0.15 

101269.7 

38150 

1100 

-0.20 

199973.1 

38150 

1100 

-0.25 

101160.3 

38150 

1100 

-0.30 

101183.6 

38150 

1100 

-0.35 

1 

101079.4 

38150 

1100 

-0.40 

101068.8 

38150 

1100 

-0.46 

101176.1 

38150 

1100 

-0.67 

101073.9 

38150 

1100 

-0.125 

200808.4 

38150 

1100 

-0.075 

201008.1 

38150 

1100 

-0.175 

101249.8 

38150 

330 

-0.175 

99781.7 

152600 

330 

-0.125 

99799.6 

152600 

Remarks 


t  (Estimated  shifted 
*  (Transfer  function 


frequency) 
changed  to 


38150) 


:  ^  (Transfer  function) 
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Table  Mean  properties  of  junction  vortex  structure  on  different  Re^ 


Junction  Vortex  Structure 


;  ”  330„ 

1100' 

6700 

_.0.283 

0.190 

6.5 

,  Vm..* 

-OAlVr^tf 

-0.48l7re/ 

~0.49C^reA 

(-0.125,0.0078) 

(-6.075,0.0032) 

(-0.153,0.0047) 

(-0.15,0.64) 

(-0.125, 0.032) 

(-0.20,0.05) 

-0.175... -6.10- 

-0.20... -r  0.075 

'  Y/tt2  ,1 

0.0...  0:0533 

0.6...  0.0634 

(backflow) 

^  {Vmuz  location) 

^  (vortex  core -location) 
*  (bimodal  flow  range) 
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Table  7.  File  list  related  to  LDV  measurements 


File  Name 

File  Size* 

Remarks 

ldvkim,c 

2,822 

main  program  for  LDV  data  reduction 

Idvavg.c 

518 

routine  for  mean  values 

Idvfinam.c 

1,134 

filename  determination 

Idvinit.c 

204 

intialization 

Idvinput.c 

820 

input  file 

non.dim.c 

204 

Idvpdf.c 

1,127 

routine  for  probability  density  function 

swfa.c 

7,231 

main  program  for  switching  frequency  analysis 

swffln.c 

344 

file  name  determination  routine 

fourl.c 

966 

data  smoothing  routine^ 

f-mem.c 

327 

free  dynamically  allocated  memory 

getonem.c 

2,001 

allocate  memory  dynamically 

hstgrm.c 

3,040 

find  a  velocity  histogram 

swfinit.c  - 

521 

initialize  the  arrays 

swfinput.c 

1,535 

input  file 

nrutil.c 

3,706 

t 

realft.c 

1,006 

t 

smooft.c 

739 

t 

smooth.c 

592 

histogram  smoothing  routine 

shJist.c 

2,986 

display  program  of  output  histogram 

switch.c 

11,466 

analysis  progreim  for  switching  frequency 

*  (in  bytes) 

^  (from  Numerical  Recipes  in  C) 


Table  8.  File  list  related  to  PIDV  motion  analysis 


File  Name 

File  Size* 

Remarks 

imageflo.c 

22,625 

main  program 

centroid.c 

3,346 

centroid  search  program 

chkerfrm.c 

692 

program  to  check  error  frame  LUT 

corresp.c 

4,453 

correspondence  checking  routine 

correspO.c 

11,631 

correspl.o 

923 

corresp2.  c 

3,383 

datinSd.c 

2,019 

input  file 

distanl.c 

684 

distance  between  bubbles  in  an  image  plane 

distan2.c 

684 

distance  between  bubbles  in  an  image  plane 

enhanim.c 

2,519 

image  enhancement  algorithm 

filename. c 

643 

a  routine  generating  file  name 

group.c 

539 

a  routine  to  check  grouping  of  bubbles 

grp.avg.c 

949 

a  routine  to  solve  for  the  average 

X  or  Y  location  in  the  same  group 

g.trans.c 

2,043 

transformation  from  2-D  to  3-D 

lubksb.c 

373 

back  substitution  routine^ 

ludcmp.c 

1,114 

LU-decomposition  routine^ 

match.c 

5,813 

matching  algorithm 

ml-coef.c 

3,803 

a  routine  for  a  corresponding  line 

m2_coef.c 

3,803 

nrutil.c 

3,497 

t 

reinit.c 

2,109 

initializing  arrays  in  the  program 

renum.c 

7,127 

bubble  elemenet  numbering  routine 

thdval.c 

485 

a  routine  for  threshold  value 

thrshold.c 

1,322 

velinst.c 

1,733 

a  routine  for  instantaneous  velocity 

veLavg.c 

1,853 

a  routine  to  retrieve  the  lost  bubble  information 

vthrd.c 

838 

a  routine  for  varying  threshold  value 

*  (in  bytes)  :  ^  (from  Numerical  Recipes  in  C) 


Table  9.  Some  statistics  from  flow  visualization  results. 


9.07cm/a,  {Re^ 

=  330) 

Fraction  of  a  backflow  mode 

0.78  ^  0.23 

Strouhal  number 

0.28  *  0.07 

Average  lapse  time  between 
the  formation  of 
the  same  events 

t 

4.26  *  0.39(scc) 

Average  duration  time  of 
a  primary  vortex  before 
the  separation  vortex  forms 

0.85  ^  0.08(5ec) 

Average  total  time  for 
"leap-frogging”  process 
to  occur 

1.87  ±  0.16(sec) 

^(estimated  dA  XfT=i  -0.125) 
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Table  10.  Comparison  of  dimensionless  frequency  groups 

obtained  by  other  studies 


present  study 

Devenport 
and  Simpson 
(1990a) 

Oilmen 
ti  al 
(1991) 

Hasan 
Casarella 
Rood  (1985) 

Re$ 

330 

1,100 

6,700 

4,450 

6,560 

SIT 

0.283 

0.190 

0.502 

0.261 

0.5 

0.068 

0.08 

0.02  (0.048)* 

0.03 

0.0032  -  0.032 

0.34 

0.40 

0.10  (0.24)* 

0.15 

0.016  -  0.16 

_ _ 

0.096 

0.076 

0.05  (0.12)* 

0.039 

0.008  -  0.08 

*  (most  probable  value) 
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FIGURES 


FIGURES 


Figure  1.  Streamline  pattenui  in  the  plane  of  symmetiy  of  small  obstacle  8q>ar 

rated  flow  showing  (a)  six  vortex  system  and  (b)  four  vortex  system. 
S,  separation  line  ;  A,  attachment  line  ;  k,  height  of  the  obstacle; 
d,  oncoming  boundary  layer  thickness  (from  Sedney  and  Kitchens 
(1975)) 
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The  water  tunnel  facility 


Figure  S(a).  Hot  film  measurements  for  the  turbulence  level(%  of  Uff)  in  the 
test  section  with  a  wing-body  junction  model  mounted  for  low 
speed  case  :  *  (looking  upstream).  (X,y,Z)  coordinate  uses 
a  right-handed  system  with  origin  at  the  comer  between  leading 
edge  of  the  wing  and  false  wall. 
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{Location*) 


{X  =  -81.3cm) 


XXX 

(48.72, 1U4)  (48.710)  (48.7a^<tlLat) 


X  X 

(20^18.24)  (2048,0) 


y  (cm) 


z  (cm) 


X  X 

(18J4,18J4)  (18J«,0) 


^  o,eie 

0048,-1B.2t) 


(lOM'lfcai) 


{X  =s  -61.0cm) 


{X  =  —30.5cm) 


0.78T  1,112- 


0,102  0,811- 


0.768  0.632- 


Figure  8(b).  Concluded  —  high  speed  case 
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Side  wall  location,  with  boundary  layer  displacement  thickness  cor¬ 
rection 


# 


(a)  Conventional  time-line  method 


(b)  Continuous  sheet  method 


Figure  7.  Examples  of  hydrogen  bubble  techniques 
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# 


(c)  Combined  time-streak  marker 


(d)  Inverted  time-line  method 


Figure  7.  Concluded 
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Figure  8.  Hydrogen  bubble  probe 


Hydrogen 


Schematic  diagram  of  hydrogen  bubble  generator  circuit 
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ptics  and  Data  Acquisition  System 
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Figure  ll.(b).LDV  setup 


Figure  IS.  Fringe  spacing  measurement  across  the  probe  volume 


(a)  Global  coordinate  system 


(238,191) 


(0.0) 


(b)  Image  plane  coordinate  system 


Figiire  14.  Coordinate  system 


FIGURES 


Mirror  Box 


6. (a).  Image  capturing  system  setup  for  Kodak  Ektapro  1000  motion  an- 


Physical 

Object 

Point 


Image 

Point 


Mathematical 

Model 

Sensor  Parameters, 
Imaging  System 
Orientation 


Figure  17. (a).  Mathematical  model  relating  3-D  object  points  to  image  points 


System 


Figure  17,  (b).  Generalized  pinhole  camera  model 


7IGUR2S 


i-2  i-1  i  i+1  i+2 
(a)  8  point  window 


i-3  i-2  i-1  i  i+1  i+2  i+3 
(b)  16  point  window 


j+2 

j+1 

j 

j-1 

j-2 


Figure  20.  Windows  for  local  averaging  technique 
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riOURBS 


.  (y  -  ymln  ) 

where  D  =  Int  ^  R  — — — - p  +  0.5 

max  ^mln  ' 

R  is  the  range  of  the  gray  level 
y  Is  the  relative  frequency 
D  Is  the  scaling  operator. 


Figure  21.  Algorithm  for  histogram  equalization  technique 
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Figure  22.  Examples  of  bilevel  images 


Identification  of  the  bubble  elements  for  centroid  search 


i-th  element 


Splitting-Merging  Detection  Algorithm 


Algorithm  for  the  detection  of  the  split-merging  bubbles 
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(d)  t  =  0.3  (sec) 


Figure  26.  Continued 


♦—  3i 

< —  d 

Cu  C2 


♦—  $1 

4—  d 

-Ci,C2 


FIGURES 


139 


d,S3 
■  32 

■CuC2 


14? 


9 


V2  Cl  V3 

(v)  t  =  2.1  (sec) 

Figure  26.  Concluded 
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' - V - '  . . . ' 

ST  P 

(b)  t  =  0.1  (sec) 

Figure  27.  Sequence  of  continuous  bubble  sheet  pictures  showing  strong  inter¬ 
action  of  two  vortices  near  the  wall 
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S  T  P 
(j)  t  =  0.9  (sec) 


Figure  27.  Continued 
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Figure  30.  Turbulent  normal  stress  distribution  across  the  boundary  layer  at 
XjT  =  —2.15  for  Ree  =  330,  Reg  =  1100  (LDV  mecisurements)  and 
Rtg  =  6700  (Hot-wire  measurements  by  Devenport  and  Simpson 
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FIGURES 


Figure  31.  Turbulent  normal  stress  distribution  across  the  bo 
X/T  =  —2.15  for  Reg  =  SSO^Reg  =  1100  (LDV  mee 
Reg  —  6700  (Hot  wire  measurements  by  Devenpor 


Estimation  of  the  local  friction  coefficient  at  X/T  =  —2.15  using 
the  “cross  plot  of  log  region”  procedure  for  Ree  =  1100 


Figure  33.  Law  of  the  wall  plot  at  XfT  =  -2.15  for  Ree  =  1100:  Ut  is  cal¬ 
culated  using  the  “cross  plot  of  log  region”  procedure  assuming 
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Figure  35.  Estimation  of  local  friction  coefficient  at  XfT  =  -2.15  using  the 
“cross  plot  of  log  region”  procedure  for  Ree  =  330 
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Y/6 


Figure  41.  Meaii  velocity  profile  comparison  for  three  Reynolds  numbers  at 
each  measurement  location  {Ree  =  330  and  Reg  =  1100  —  present 
data  :  Reg  =  6700  —  Devenport  and  Simpson  {1990b)) 
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U/Uref 


FIGURES 


Vu2/Uref 


10-4  2  10-^  ^  10-2  ^  10-'  ^  10°  ^ 


Y/6 

Figure  42.  U-component  turbulent  fluctuations  for  three  Reynolds  numbers  at 
each  measurement  location  [Ree  =  330  and  Ree  =  1100  —  present 
data  :  Re$  =  6700  —  Devenport  and  Simpson  (1990b)) 
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Figure  42.  Continued 
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Contour  plot  for  the  upstream  of  the  nose  of  the  wing-body 
junction  for  Re$  =  330 


Contour  plot  for  the  upstream  of  the  nose  of  the  wing-body 
junction  for  Ree  =  1100 


(a)  for  Reg  =  1100  at  X/T  =  -0.10 


Figure  44. 


Velocity  histograms  across  the  boundary  layer 
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Figure  44.  Continued 
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(c)  for  Ree  =  1100  at  X/T  =  -0.150 


Figure  44.  Continued 
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(d)  for  Ree  =  1100  at  X/T  =  -0.175 
Figure  44.  Continued 
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(e)  for  Ree  =  1100  at  XjT  =  — 


Figure  44.  Continued 
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(h)  for  Ree  =  330  at  XfT  =  -O.li 


Figure  44.  Continued 
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(i)  for  Ree  =  330  at  X/T  =  -0.175 
Figure  44.  Continued 
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(j)  for  Reg  =  330  at  X/T  =  -0.: 
Figure  44.  Concluded 
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Figure  46,  Continued 
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Figure  46.  Concluded 
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Figure  47.  Flatness  factor  plots  across  the  shear  layer  for  Re$  =  330 
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Figure  47.  Continued 


Figure  48.  Skewness  factor  plots  across  the  shear  layer  for  Ree  =  330 
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Figure  48.  Continued 
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Figure  50.  FIR  low  pass  filter  design 
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Figure  50.  Continued 
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FIR  Low  Pass  Filter  Design 
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Figure  50.  Continued 
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Figure  50.  Concluded 
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Figure  64.  Continued 
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Variations  ofjthe  average  switching  time  {Tau>)  in  the  bimodal  zone 

for  each  flow  condition 
« 


the  average  switching  frequency  {fsw)  in  the  bimodal 
flow  condition 


Contour  plots  of  the  estimated  fraction  of  backflow 
the  bimodal  zone  for  each  flow  condition 
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Figure  59.  Continued  :  (m)  t  =  to  +  0.096(5ec),  (n)  t  =  to  +  0.104(c>ec),  (o) 
t  =  to  +  0.112(5ec),  and  (p)  t  =  to  +  0.l20(i:ec) 
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Figure  59.  Continued  :  (q)  i  =  to  +  0.12S(5ec),  (r)  t  =  t 
t  =  to  +  0.144(5ec),  and  (t)  t  =  ig  +  0.152(5ec) 


+  0.136(5ec),  (s) 
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FIGURES 


Figure  59.  Concluded  :  (u)  t  =  to  +  0.160(5ec),  (v)  t  —  to  +  0.16S(^'cc),  (w) 
t  =  to  +  0.176(5ec),  and  (x)  t  =  to  +  0.1S4(o>cc) 
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Hydrogen  Bubble  Wires 
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3-D  Instantaneous  Velocity 
Vector  Field  around  a  Wing 


Figure  60. 


Sequence  of  instantaneous  three-dimensional  velocity  vector  field 
in  the  nose  region  of  the  wing-body  junction  (meastired  by  PIDV) 
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figures 


(a)  0-th  bubble  element 

Figure  61.  Three-dimensional  instantaneous  velocity  components  and  loca¬ 
tions  of  the  individual  bubble  elements  generated  from  the  hydro¬ 
gen  bubble  wire  ;  the  results  obtained  by  tracing  the  hydrogen 
bubble  patches  generated  at  the  first  wire  are  provided  from  the 
bottommost  element  to  the  top. 
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(b)  i-st  bubble  element 
Figure  61.  Continued 
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(c)  8-nd  bubble  element 


Figure  61.  Continued 
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(e)  4-th  bubble  element 
Figure  61.  Continued 
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(f)  5-th  bubble  element 
Figure  61.  Continued 
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(h)  7-th  bubble  element 
Figure  61.  Continued 


FIGURES 


(a)  Primary  vortex  (P.V.)  is  formed 
upstream  of  the  wing  (20%). 


(b)  Secondary  separation  vortex  (S.V.) 
with  the  same  direction  as  the  P.V. 
developed  upstream.  Another 
counter-rotating  vortex  exists 
between  two  strong  vortices  (15%). 


Fluids  in  low  shear  region 
separates  from  the  wall. 
Lifted-up  fluids  create 
rotational  disturbances  in  the 
outer  layer  (15%). 


S.V.  moves  downstream  while 
the  disturbances  created  upstream 
begin  to  roll  up  over  a  large 
vortical  structure.  Leap-frogging 
of  multiple  vortices  adds  more 
strength  to  the  primary  vortical 
structure  (21%). 


(e)  While  the  strengthened  P.V. 
is  stretched  around  a  wing, 
accelerated  flow  is  induced 
in  this  region.  The  inrush 
of  high  momentum  fluids  energizes 
this  region.  Flow  is  stabilized 
at  this  instant  (8%). 

Figure  62.  Descriptive  model  for  the  sequence  of  flow  events  in  the  nose  re¬ 
gion  of  a  wing-body  junction.  The  percentage  in  the  parenthesis 
represents  the  approximate  time  proportion  of  the  event. 
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unstable  flow 


(f)  After  a  while,  flow  in  this  region 
becomes  very  unstable  (7%). 


(g)  Some  time  later,  a  weak  large 
vortical  structure  is  formed  again 
upstream  of  the  wing. 

A  cycle  of  (a)  to  (g)  is 
repeated  building  up  time- 
dependent  characteristics  of 
a  junction  flow  (14%). 


Figure  62. 


Concluded 


Variations  of  Strouhal  numbers  (Sir)  in  the  bimodal  zone  for  each 
flow  condition 
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APPENDIX  A.  UNCERTAINTY  ANALYSIS 


Suppose  that  x,y,. . .  ,z  are  measured  with  uncertainties  Sx,  Sy,. . . ,  6z,  and  the 
measured  values  used  to  compute  the  function  q(x,  y,. ..  ,z).  If  the  uncertainties  in 
a;,  y, . . . ,  z  are  independent  and  random,  then  the  uncertainty  in  q  is 

=  +  +  (^.1) 

In  any  case,  it  is  never  larger  than  the  ordinary  sum 

p-  6y  +  ---+  —  Sz  (A.2) 

ox  by  bz 

In  case  of  the  random  errors,  wc  can  repeat  our  measurements  in  order  to 
achieve  reliability  of  the  experiments,  and  the  measured  data  set  now  can  be  ana¬ 
lyzed  statistically. 

For  the  random  variable  x  with  N  independent  measurements  {xi,X2, . . .  ,X!^), 
the  best  estimate  would  be  the  average. 

Xbest  =  X  (A. 3) 


Actual  imcertainty  in  random  variable  x  with  N  measurements,  however,  is  using 
another  quantity,  which  is  the  standard  deviation  of  the  mean  that  is  defined  by 

Thus,  the  value  of  the  random  variable  x  can  be  expressed  cis 

(value  of  r)  =  X  ±  crj-  (-^T) 


Here,  the  value  of  x  expressed  in  Equation  (A. 7)  has  the  confidence  level  of  68%, 
which  means  that  only  68%  of  any  subsequent  measurements  of  x  are  expected  to 
fall  in  the  range  x  ±  crj.  If  the  final  value  is  written  by 

(value  of  x)  =  X  ±  2aj  (^-8) 

or 

(value  of  x)  =  X  ±  ScTj,  (^-9) 


the  confidence  level  of  each  value  would  be  95%  and  99%. 

The  present  study  used  confidence  level  of  68%  to  compute  the  uncertainties  of 
the  individual  measured  variables,  and  then  Equation  (A.l)  was  applied  to  obtain 
the  combined  uncertainties  of  the  target  quantity  q. 
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APPENDIX  B.  Q  -  R  ALGORITHM 


Suppose  we  are  given  an  overdetermined  linear  equation  of  the  form 


Ax  =  b  (B.l) 

where  6  is  m  x  1,  s  is  n  x  1,  and  A  is  an  m  x  n  matrix  of  rank  m. 

Q  -  R  method  is  based  on  the  decomposition  of  the  m  xn  matrix  A  into 

A=^QR  (5.2) 

where  Q  is  m  x  m.  and  satisfies 

Q'^Q  =  D  (5.3) 

D  is  the  diagonal  matrix  with  non-zero  diagonal  elements.  R  is  upper  triangular 
n  X  n  matrix  with  diagonal  entries  rkk  =  1.  (Notice,  therefore,  R~^  exists.)  Given 
Equation  (B.l)  and  Equation  (B,2),  we  may  reformulate  the  least-squares  solution 
as 

R^Q'^{b-Ax)  =  0.  (5.4) 

Since 

Q'^A  =  Q'^QR  =  DR  (5.5) 

the  preceding  equation  may  be  rewritten  as 

R'^{Q'^b-DRx)  =  0.  (5.6) 

Thus,  we  arrive  at 

DRx  =  Q^b 

or 

Rx  =  D-^Q^b.  {B.l) 

Thus,  if  this  equation  were  formulated,  we  could  solve  for  x  via  back  substitution. 
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APPENDIX  C.  m,’  COEFFICIENTS 


The  coefficients  m,’s  in  Equation  (5.15)  can  be  derived  by  arithmetically  ex¬ 
panding  Equation  (5.14)  as  follows, 

+  ai32<^  +  O’lizP  +  0134’^^ 

m2  =  —  <2112^^  —  0113-^  ~  O114T 

^3=  o.\z\E  o.\z2K o.\zzQ 
Tn.x  =  — flinA  —  0112^’  —  aiisAf  —  aii4>S 
W5  =  013iF  -b  ai32-b  +  0133-R  +  0134-^ 

mg  =  — ajjiD  —  aii2*I  —  onsF  —  o.\\^V 
m-j  =  — QiiiC  —  0112/  ~  o,\\zO  —  a\\\U 

niQ  =  — oiiiF  —  0112/1  —  O113Q  —  diiA^V 

mg  =  ~aui/'  ~  0112/'  —  a,\\zR  —  clwaX. 

Here,  the  coefficients  (mi ,  •  •  • ,  mg)  are  derived  for  the  coefficients  of  the  line  of 
correspondence  in  the  upstream  view  image  plane  for  a  point  in  the  side  view  image 
plane.  The  constants  A  thru  X  are  functions  of  camera  calibration  coefficients  and 
expressed  in  the  following  pages. 
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